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ABSTRACT 

Under  this  AFOSR  GameChanger  we  carried  out  groundbreaking  research  by  developing 
methods  and  materials  to  realize  small  and  conformal  structurally  integrated  antennas  for 
unmanned  airborne  vehicles  (UAVs).  We  are  proud  to  report  several  groundbreaking 
developments  in  all  areas.  Among  them  we  (a)  developed  a  conformal,  load-bearing  antenna 
technology  by  embroidering  e-fibers  on  polymer  composites,  (b)  Characterized  performance  of 
the  e-fibers  in  terms  of  core  and  specific  strength  for  load-bearing  applications,  (c)  Developed  a 
new  principles  DFT  method  for  systematic  search  of  ferromagnetic,  low  conductivity  nanowires, 
(d)  Developed  a  strategy  for  increasing  permittivity  and  permeability  proportionately  using 
magnetic  nanowires  combined  with  conventional  high  contrast  materials  for  antenna 
miniaturization,  (e)  demonstrated  theoretical  feasibility  of  laminate  EM  materials  that  alter  wave 
propagation  characteristics,  critical  to  future  miniaturization  efforts,  (f)  Developed  a  multi-time 
scaling  model  and  algorithm  for  coupling  mechanical  and  EM  solutions  in  real  time,  critical  to 
multiphysics  coupled  EM-Structural  analysis,  (g)  Extended  our  new  domain  decomposition 
algorithm  to  3-D  space  and  characterized  more  complex/intricate  details  on  large  platforms  such 
as  e-fibers  embroidered  on  UAV  platforms.  In  addition  to  this  report,  our  ground  breaking  work 
under  this  program  has  been  published  in  over  120  journal  papers  and  over  70  conference 
summaries. 


1 


SUMMARY  OF  ACCOMPLISHMENTS . 3 

Antennas  and  RF  Systems: . 3 

Materials  and  Composites: . 3 

Multiphysics/Multiscale  Analysis/Design/Optimization: . 3 

Integration  and  Transition: . 3 

PROF.  VOLAKIS’  GROUP  SUMMARY: . 4 

E-Fiber  Embroidered  on  Polymer  Composites . 9 

PROF.  ROJAS’  GROUP  SUMMARY . 18 

PROF.  KOTOV’S  GROUP  SUMMARY . 32 

Research  Summary . 41 

PROF.  SHTEIN’S  GROUP  SUMMARY . 42 

PROF.  GHOSH’S  GROUP  SUMMARY . 44 

PROF.  BECHTEL  ’S  GROUP  SUMMARY . 54 

PROF.  CARMAN’S  GROUP  SUMMARY . 61 

Research  Summary . 62 

Analytical  Wave  Propagation  Modeling: . 63 

Development  of  effective  laminate  properties: . 68 

Fabrication,  Characterization  and  Testing  of  ME  materials . 70 

PROF.  JAMES’  GROUP  SUMMARY . 77 


2 


Summary  of  Accomplishments 

Antennas  and  RF  Systems: 

•  Developed  embroidered  E-fiber  -  polymers  for  conformal,  load  bearing  applications 

•  Developed  ,  built  and  characterized  3-D  volumetric  E-fiber  spiral  antenna  in 
collaboration  with  U  of  Michigan. 

•  Extended  the  2-D  Domain  Decomposition  Algorithms  to  3-D  domain  and  demonstrated 
their  convergence  for  nanowires  embedded  within  composites. 

Materials  and  Composites: 

•  Optimized  process  for  e-fiber  fabrication 

•  Characterized  performance  of  e-fibers  in  regards  core  material  and  specific  strength, 
conductance,  twist  multiplier  etc. 

•  Developed  conductive  textiles  based  on  CNT  and  copper  sputtering  on  Kevlar 

•  Fabricated  thin  film  magneto-electric  structures 

•  Shown  theoretical  feasibility  of  laminate  ME  materials  to  significantly  alter  wave 
propagation  characteristics 

•  Developed  numerical  techniques  to  characterize  magnetic  semiconductors  and 
electromagnetic  materials. 

•  Developed  strategy  for  increasing  sr  pr  proportionately  using  magnetic  nanowires 
combined  with  conventional  high  dielectric  constant  materials 

•  Developed  low  loss  design  principles  by  prediction  of  breakdown  of  the  single  domain 
state  in  nanowires 

•  Developed  a  new  first  principles  DFT  method  for  systematic  search  for  ferromagnetic, 
low  conductivity  nanostructures,  based  on  the  concept  of  objective  structures 

Multiphysics/Multiscale  Analysis/Design/Optimization: 

•  Developed  Coupled  Electro-Magnetic  and  Dynamic  Response  codes 

•  Developed  a  multi-time  scaling  model  and  algorithm  for  coupling  mechanical  and  EM 
solutions  in  real  time 

•  Introduced  high-performance  parallel  computational  framework 

•  Multi-Scale-Super-Element  which  exploits  shape  functions  derived  from  semi-analytical 
solutions  in  improving  accuracy  and  enabling  parameter  studies  for  design  and 
optimization  problems 

Integration  and  Transition: 

•  Develop  prediction  formulae  for  dielectric  constants  as  a  function  of  temperature  and 
stress  (Correlated  these  parameters  to  antenna  resonance  and  predicted  resonance  shift) 

•  Developed  model  for  antenna  performance  under  various  conditions 
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Research  Summary 

A  key  development  under  the  Gamechanger  program  is  a  new  breakthrough  technology  based  on 
embroidering  E-fibers  on  polymer  substrates  for  conformal  and  loadbearing  antennas.  As  is 
already  known,  polymers  provide  for  conformal,  flexible,  light-weight  substrates  that  can  be 
readily  reinforced  for  load-bearing  antennas  on  small  UAVs  and  body  worn  applications. 
Electrically  conductive  fibers,  comprised  of  thin  metal  coated  polymer  strands  (Zylon,  Kevlar 
etc.),  are  extremely  flexible.  They  are  also  easy  to  stitch  onto  polymer  substrates  forming  an 
embroidering  structure.  Also,  the  strands  are  extremely  conductive  (comparable  with  traditional 
conductors).  Therefore,  there  is  no  noticeable  compromise  on  the  antenna  performance  and 
efficiency. 

In  this  technology,  polymer-ceramic 
composites  serve  as  the  host  media 
(see  Figure  1).  E-fibers  are  then 
embroidered  on  the  polymer  to 
realize  planar  and  volumetric 
antennas.  Similarly,  these  fibers  can 
also  be  used  to  form  a  grid  ground 
plane  for  structural  integration  with 
the  platform.  In  other  words,  these 
conductive  and  extremely  strong 
fibers  can  be  stitched  into  the 
platform  for  higher  structural 
stability.  An  optional  layer  of  non- 
conductive  fiber  such  as  Kevlar  or 
Zylon  can  also  be  embedded  inside 
the  polymer  composite  to  enhance 
structural  ruggedness  of  the  overall 
antenna  structure  (see  Figure  1). 

A  key  component  of  this  technology  is  its  low-cost,  scalable  and  automated  fabrication  process 
for  embroidering  e-fibers  on  polymer  composites.  This  is  a  game-changing  technology  whose 
realization  will  provide  a  new  paradigm  in  load-bearing  electronics  and  RF  front  ends.  At  this 
moment,  we  have  only  demonstrated  the  potential  of  the  E-fiber  embroidering  technology.  The 
next  obvious  and  critical  step  is  to  control  printing  precision  and  realize  multilayer  low-loss 
(high  efficiency)  front  end  antenna  arrays  with  electronic  back-ends. 

The  unique  aspects  of  our  e-fiber  embroidering  technology  are: 

•  Straightforward  and  easy  to  scale  fabrication  in  room  temperature 

•  Flexible,  light-weight  conformal,  load  bearing  and  extremely  durable  structure 

•  Chemically  configurable  polymer  composite  for  higher  permittivity  and  permeability 
metamaterials  for  antenna  miniaturization 


loadBearlna  Polymer  Based  Antenna  Concept 


(oorfmttv*  to#  A/it*nrw  PflMtatf 

(UMHlutt**.  f !•«**>,  umfomial  bourn  a 
tout**.  CN1  or  •  f#«d 


Cw«nk  CorapotM*) 
IW«Ot,  UN*  IOM 


Rrtntutcmn*  toiwtucthr* Ground  Win* 1  i  /  // 
(<«nb4>  wmvnI  toto  ft**  %Un  ot  tho  Unw  tuie  | 

Conduct**  ttrondt,  Amb*n1r«nd'.  UM*h*~f»b*r» 


Figure  1  Embroidered  E-fiber  polymer  composite  for  load-bearing 
conformal  antennas 
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•  Printed  E-fiber  can  be  flexed  without  concern  about  detachment  and  performance 
deterioration  due  to  vibration  and  stretching. 


In  the  final  year  of  the 
GameChanger  program,  we 
continued  to  work  closely  with 
Max  Shtein  of  University  of 
Michigan  to  develop  volumetric 
antennas  using  copper  coated 
Kevlar  yarns  stitched  into 
polymer  substrate  (see  Error! 
Reference  source  not  found.). 
Our  previous  efforts  in  coating 
Kevlar  with  copper  resulted  in  a 
highly  conductive  yarn  weaved 
into  a  polymer  substrate.  The 
final  weaved  antenna  is 
extremely  flexible  and 
conformal  to  any  surface  of 
interest. 


Figure  2  qiiw1  »w*eiiwai  printed  by  U  of 

Michigan  (a)  sample  of  encapsulated  conductive  yarn,  (b) 
patterned  PDMS  substrate,  (c)  example  of  w  eaving  and 
gluing  (inset)  (d)  completed  flexible, woven  antenna. 


Error!  Reference  source  not 
found.  displays  the 

performance  of  a  stitched 
volumetric  spiral  antenna  with 
that  of  a  traditional  one  made  of 
copper  wires  one  a  ceramic 

substrate.  We  notice  the 
remarkably  good  performance 
at  frequencies  above  700  MHz. 

At  lower  frequencies,  the 
weaved  e-fiber  antenna 
experiences  lower  gain, 

primarily  due  to  resistive  losses 
on  the  copper  coated  Kevlar.  In 
other  words,  at  lower 

frequencies,  more  actual  length 
of  the  wire  is  active  to  achieve 

radiation,  leading  to  higher  resistive  losses.  Such  losses  become  less  pronounced  at  higher 
frequencies  since  the  antenna  resonance  is  achieved  using  relatively  shorter  e-fiber  lengths. 
Clearly  higher  conductivity  is  needed  for  the  copper  coated  Kevlar  for  low  frequency 
applications. 


Vofctnwlrtc  E4lb«f  *r«l*nrv» 

*  PDMS  Rfc  T  ♦  dided  ♦ 

£  Abe* 


Figure  3  RF  performance  of  the  sdtched  volumetric 
spiral  antenna  compared  with  that  of  metallic  ceramic 
antenna 


In  a  parallel  effort,  we  made  remarkable  progress  in  domain  decomposition  methods  for  fast  and 
convergent  computations  associated  with  domains  involving  antenna  details,  particularly  feeds  or 
small  wire  and  surface  features.  Last  year,  we  demonstrated  our  novel  fully  overlapping  domain 
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decomposition  approach  based  on  2-dimensional  (2D)  geometries.  Using  radiation  and  scattering 
examples,  we  illustrated  its  superior  accuracy  and  fast  convergence.  Particularly,  our  overlapping 
DDM  was  found  suitable  for  problems  involving  multi-scale  modeling  of  composite  structures 
and  embedded  antennas.  It  achieves  fast  convergence  by  decoupling  the  ultra-dense  mesh  for 
small  features  from  the  homogeneous  background  mesh.  In  the  last  year  of  the  program,  we  are 
happy  to  report  the  extension  of  the  overlapping  DDM  to  3 -dimensional  (3D)  problems. 
Technical  details  on  this  achievement  are  provided  in  the  subsequent  section. 


Embroidered  E-fiber  Technology 


E-Fiber  Embroidering  on  Polymer  Composites 

To  demonstrate  the  embroidering  concept,  we  used  Amberstrand®  silver-coated  Zylon  fibers  as 
the  E-fiber  conductor  (see  Figure  3).  This  e-fiber  had  a  resistance  of  0.8  Q/m  in  addition  to 
superior  mechanical  strength  and  flexibility  due  to  the  core  Zylon  fiber.  As  for  the  RF  substrate, 
we  employed  polydimethylsiloxane  (PDMS)  mixed  with  rare-earth  titanate  ceramics  (RET).  The 
PDMS-RET  composite  was  comprised  of  90  vol%  of  PDMS  and  10  vol%  of  RET.  With  this 
composite,  we  achieved  a  dielectric  constant  of  sr=4.15  and  a  low  loss  tangent  of  tan  S  <0.01. 


Bottom' 

1st  reinforcing  fabric 


We  developed  a  fabrication  process 
to  integrate  the  e-fibers  onto  the  polymer 
substrate  (see  Figure  3).  The  substrate  is  also 
reinforced  by  non-conductive  fabric  layers  to 
improve  its  load-bearing  capability  as  shown 
in  Figure  1.  By  employing  the  process  in 
Figure  3,  we  fabricated  a  sample  patch 
antenna  and  the  details  of  the  fabrication 
process  will  be  presented  at  the  conference. 

Although  the  E-fiber  patch  does  not  have  a 
physically  continuous  conductive  surface,  it 
can  be  regarded  as  an  electrically  continuous 
surface  since  it  was  tightly  knitted.  The  resolution  of  the  e-fiber  mesh  was  much  less  than  1/20  of 
the  wavelength  of  operation. 


Medium  W  Top 

PDMS/ceramic  composite  2nd  reinforcing  fabric 


Add  PDMS/ceramic 
to  Back  side 


I  Stitch  ,  /Ambers^and® 

mm  w  e-fiber 


Substrate~5mm  thick  Substrat^mm  thick 

Figure  3  E-fiber  Polymer  Substrate  Fabrication  process 


As  displayed  in  Figure 
1,  an  additional  e-fiber  surface 
can  also  be  used  to  form  the 
ground  plane.  Further,  non- 
conductive  fibers  such  as 
Kevlar,  Zylon  can  be  employed 
in  the  middle  layer  to 
strengthen  the  substrate  and 
allow  for  load  bearing 
capability. 


E- fiber  Automated 
Embroidery  Process 


•  Electronically  conductive  fibers 

-  Ag  coated  Zylon  fibers,  332  threads  in  one  strand. 

-  Thick,  easily  spread:  not  good  for  direct  sewing. 


•  Embroidered  E-fiber  array 

-  E-fibers  do  not  go  through  the  needle  or  substrate. 

-  Embroider  the  substrate  with  the  antenna  side  down, 
use  assistant  yams  to  "couch"  down  E-fibers. 

-  Digitize  the  route  for  yams  to  obtain  the  designed 
antenna  pattern. 


Yam  stitching  route 


◄ 


Assistant  yams 


Design 


Embroidered  antenna 


E-fibers 

(one  continuous  strand  of  16.6m!) 


Antenna  side 


Bobbin 


E-fibers-> 


Figure  2  E-fiber  Embroidery  Process 


We  closely  worked  with  embroidery  experts  to  identify  the  best  feasible  stitching  technique  for 
the  type  of  e-fiber  we  have.  It’s  important  to  note  that  the  stitching  yam  used  for  sewing  and 
embroidering  are  specially  made  by  twisting  the  fiber  into  more  robust  form  to  be  able  to 
accommodate  large  number  of  frictions  during  the  process.  Given  that  our  e-fiber  is  a  silver 
coated  fiber,  any  friction  during  the  stitching  process  will  degrade  RF  performance  by  reducing 
the  fiber’s  conductivity.  We  therefore,  identified  the  bobbin  stitching  technique  as  the  ideal 
approach  for  embroidering  e-fiber  on  the  polymer  substrate. 

The  above  embroidery  process  is  very  scalable  and  computer  automated  as  shown  in 
Figure  2.  The  process  starts  with  exporting  the  antenna  design  from  widely  used  RF  software 
tools  such  as  HFSS,  FEKO  etc.  to  graphics  form.  This  is  next  fed  into  a  digitization  software 
used  for  embroidering.  This  software  identifies  the  stitches  needed  to  make  this  pattern  as 
accurate  as  possible.  The  digitations  file  is  next  imported  to  a  computerized  embroidery  or 
sewing  machine  for  final  printing.  This  resulted  in  a  pattern  printed  on  a  thin  layer  of  polymer  of 
textile  for  further  integration  into  polymer  with  the  process  shown  in  Figure  3. 

E-fiber  RF  Patterns  Fabrication  Process 

As  shown  in  Figure  6,  the  embroidered  RF  patterns  are  integrated  on  polymer  substrate  during 
the  fabrication  process.  It  starts  with  the  mixing  of  PDMS  composite  in  room  temperature  , 
followed  by  partially  curing  process  of  the  polymer  mixture  under  elevated  temperature.  The 
mixture  is  next  tapecasted  to  form  the  substrate  with  smooth  surface.  Meanwhile,  the  backside  of 
the  embroidered  E-fiber  textile  pattern  was  cleaned  up  for  the  next  printing  process.  As  depicted 
in  Fig.  6,  the  E-fiber  pattern  was  carefully  flattened  on  top  of  the  partially  cured  substrate, 
followed  by  a  fully  curing.  We  fabricated  a  50Q.  polymer-based  embroidered  microstrip  line. 
Final  configuration  of  this  microstrip  was  highly  flexible  with  strong  adhesion  between  the  E- 
fiber  textile  and  the  PDMS.  The  entire  fabrication  process  is  low  cost  and  straight-forward, 
making  it  a  versatile  printing  method  for  various  planar  circuit  and  antenna  specification. 


lap<  ca«tir>£  polymci 


\%*emblc  inthouleKd 


Printed  embroidered 
RF  circuit _ 


Figure  6  E-fiber  RF  patterns  Fabrication  Process 
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E-fiber  RF  Performance  Evaluation 

We  next  looked  into  evaluating  performance  of  embroidered  transmission  line  and  antennas  on 
polymer  substrates. 


A.  E-Fiber  Transmission  Line  Performance 


The  RF  performance  of  the  embroidered  E-fibers  was  evaluated  using  the  50Q  transmission  line 
(TL)  structure.  As  shown  in  Fig.  7,  three  samples  were  fabricated  on  PDMS  substrates.  They  are 
composted  of  (A)  E-fiber  TL  and  copper  ground  plane  (GP),  (B)  E-fiber  TL  and  E-fiber  ground 
plane,  and  (C)  Copper  TL  and  copper  ground  plane.  Their  S -parameters  were  measured  using  an 
Agilent  N5230A  network  analyzer.  As  shown  in  Fig.  8,  while  sample  (A)  and  (C)  had  the  same 
copper  ground  plane,  the  insertion  loss  (S21)  of  the  sample  (A)  composed  of  E-fiber  TL  was 
<  0.18dB/cm  up  to  4  GHz,  only  0.04dB/cm  higher  than  that  of  the  sample  (C)  of  copper  TL.  In 
addition,  sample  (A)  exhibited  good  impedance  match  up  to  6GHz.  This  result  demonstrated  the 
high  conductivity  of  the  embroidered  E-fibers  and  low  loss  property  of  the  flexible  PDMS 
substrate  at  RF  frequencies. 

Meanwhile,  sample  (B)  showed  a  low 
insertion  loss  <  0.21dB/cm  up  to  4GHz, 
about  0.07dB/cm  higher  than  that  of  the 
sample  (C).  It  also  exhibited  a  good 
impedance  match  up  to  5GHz.  However, 
compared  with  sample  (A),  sample  (B) 
had  a  higher  IS21I  which  indicates  a 
slightly  increased  loss.  This  is  possibly 
attributed  to  the  power  dispersion  due  to 
the  physical  discontinuity  of  the 


(A)  (B)  <0 

Figure  7  Three  transmission  line  samples  (A)  E-fiber  TL 
and  E-fiber  ground  plane,  (B)  E-fiber  TL  and  copper 
ground  plane,  and  (C)  Copper  TL  and  copper  ground 
plane. 


embroidered  E-fiber  ground  plane. 


B.  Planar  E-Fiber  Patch  Antenna  Performance 


Figure  9  Simulation  model  and 
fabricated  E-fiber  patch  antenna 


PDMS-RET  substrate 


Probe  feeding 
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Referring  to  Figure  9,  we  fabricated  a  31mm  x  31mm  e-fiber  patch  on  a  56mm  x  56mm  PDMS 
substrate.  The  polymer  substrate  had  a  dielectric  constant  of  sr=3.8  and  loss  tangent  of 
tan8=0.015  at  2.2  GHz  (resonance  frequency).  As  shown  in  Figure  10,  the  measured  gain  was  5.7 
dB,  viz.  only  0,2  dB  lower  than  the  simulated  value  corresponding  to  a  perfect  electric  conductor 

(PEC)  patch.  This  clearly  demonstrates  the  remarkable  performance  of  the  e-fiber  conductors. 
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Figure  10  (a)  Return  loss  (Sn) 


(b)  Realized  gain 


(c)  Radiation  pattern 


C.  E-Fiber  Patch  on  a  Cylindrical  Surface 

To  demonstrate  performance  the  e-fiber  patch’s  on  a  conformal  surface,  the  patch  was  mounted 
on  a  cylinder  surface  (see  Figure  11).  Referring  to  the  Figure  for  the  measured  return  loss  (Sn), 
the  gain  and  the  radiation  pattern  of  the  bent  E-fiber  patch,  we  observe  excellent  agreement 
between  the  measurements  and  simulations.  As  expected,  the  curved  patch  has  lower  gain  and 
this  is  predicted  by  the  measurements  and  simulations.  Measurements  and  simulations  have  a 
mere  ldB  discrepancy  in  the  gain  as  this  is  related  to  complexity  of  this  setup  for  the  simulations. 


Frequency  (GHz) 
(b) 


Figure  8  Measured  S-parameters  of  three  transmission  line  samples:  (a)  measured  S2i  (b)  measured  Sn. 
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Figure  11  Fabricated  e-fiber  patch  antenna  mounted  on  cylinder  and  the  return  loss 


Figure  12  Measurement  RF  Performance  of  a  Cylinder  Mounted  e-fiber  patch  antenna 


I).  Embroidered  E-Fiber  Patch  Array 

A  4x1  E-fiber  array  antenna  was  fabricated,  measured,  and  compared  with  the  measured  copper 
antenna  and  simulated  PEC  antenna  of  the  same  design  (see  Fig.  13).  As  shown  in  Fig.  14,  the  E- 
fiber  antenna  had  a  resonance  frequency  of  2.31GHz  and  the  realized  gain  of  7.0  dB.  This 
resonance  frequency  matched  very  well  with  that  of  the  copper  and  PEC  antennas,  and  the 
realized  gain  was  only  ldB  lower.  On  the  other  hand,  the  discrepancy  in  the  resonance  frequency 
between  the  three  arrays  was  possibly  introduced  by  the  discontinuous  conductive  textile 
structure  that  leads  to  a  slightly  frequency  shift. 

Next,  the  performance  of  the  conformal  E-fiber  array  antenna  on  the  curvilinear  surface  were 
measured  similar  to  that  of  the  conformal  E-fiber  patch  antenna.  For  comparison,  a  Cu  array 
antenna  was  fabricated  by  prefixing  it  on  the  conformal  substrate  (see  Fig.  15).  The  measurement 
results  of  E-fiber  patch  array  were  in  good  agreement  with  both  copper  array  and  simulation,  as 
shown  in  Fig.  16.  The  realized  gain  of  E-fiber  array  was  4.6dB,  ldB  lower  than  that  of  the  Cu 
array.  In  addition,  the  increased  spacing  among  the  patch  elements  during  bending  led  to  the 
difference  between  the  radiation  patterns  of  the  conformal  and  flat  array  antennas. 

The  remarkable  results  of  the  E-fiber  transmission  line,  patch  antenna,  and  array  antennas 
confirmed  the  feasibility  of  the  proposed  embroidered  E-fiber  technology  for  emerging 
conformal  and  flexible  RF  electronics.  Moreover,  the  entire  fabrication  process  was  time  and 
cost  efficient,  and  has  a  great  scale-up  potential  for  industrial  fabrications  and  applications. 
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Figure  14  Measured  RF  performance  of  the  planar  E-fiber  antenna  array:  (a)  Sn;  (b)  realized  gain;  (c) 
H-plane  radiation  pattern 


(A)  (B)  (C) 


Figure  4  E-fiber  array  antenna  samples  mounted  on  curvilinear 
surface  (A)  HFSS  simulation,  (B)  E-fiber  array,  and  (C)  Copper 
array 
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Figure  5  Measured  RF  performance  of  the  E-fiber  antenna  array  mounted  on  curvilinear 
surface:  (a)  Sn;  (b)  realized  gain;  (c)  H-plane  radiation  pattern 


Fully  Overlapping  Domain  Decomposition  (DDM)  for  Modeling  Intricate 
Details  within  Large  Media 

In  previous  years,  we  proposed  a  new  fully  overlapping  domain  decomposition  algorithm 
for  the  finite  element  modeling  of  intricate  details  embedded  within  large  non-metallic  media. 
This  final  year,  we  further  demonstrated  the  validity  and  efficiency  of  this  novel  approach  using 
/?-refinement  convergence.  Also,  higher  order  basis  functions  and  higher  order  transmission 
conditions  were  shown  to  greatly  improve  the  accuracy  and  flexibility  of  our  method. 

The  /?-refinement  and  DDM  iterative  convergence  of  the  proposed  fully  overlapping 
domain  decomposition  method  are  compared  with  those  of  the  conventional  non-overlapping 
scheme,  using  a  canonical  PEC  spherical  scattering  problem.  In  Figure  17,  the  scattering  sphere, 
with  radius  R,  is  small  compared  to  wavelength  (R  =  / 1/8).  It  is  shown  in  Figure  17(a)  that 
fully  overlapping  DDM  approach  converges  with  respect  to  refined  mesh  density  in  almost 
exactly  the  same  rate  as  non-overlapping  scheme  and  direct  finite  element  method.  However, 
much  faster  iterative  convergence  and  much  less  iterative  time  can  be  achieved  using  the  fully 
overlapping  DDM  method  as  shown  in  Figure  17(b). 


15 


;g  -10 


I 

rflr  Vv 
Ya 

. 

-  Standard  FEM 

Non-Overlapping  DDM  (Conformal) 
Non-Overlapping  DDM  (Non-Conformal) 

1  Chimera  DDM 

■Fully  Overlapping  DDM  (Non-Conformal) 


5  10  15  20 

Number  of  Edges  per  Wavelength 


(a)  (b) 

Figure  17:  (a)  h-refinement  convergence;  (b)  iterative  convergence;  for  a  PEC  spherical  scattering  problem 

(/?  =  / 1/8). 


Up  till  now,  we  have  presented  the  properties  of  the  fully  overlapping  domain 
decomposition  approach  using  first  order  Nedelec  edge  elements  and  first  order  Robin 
transmission  conditions.  Similar  to  direct  finite  element  method  or  non-overlapping  DDM, 
higher  precision  can  be  achieved  by  using  higher  order  edge  elements  and  transmission 
conditions.  Furthermore,  the  detail  domain  boundary  can  be  placed  very  close  ( r{  «A/4)  or 
even  touching  the  antennas  when  using  higher  order  transmission  conditions.  For  instance,  in 
antenna  in-situ  design  and  optimizations,  we  frequently  need  to  place  the  antenna  on  large  non- 
metallic  platforms.  In  these  scenarios,  the  detail  domain  boundary  would  be  touching  the  antenna 
itself.  Next,  we  demonstrate  the  accuracy  and  iterative  convergence  of  an  antenna  in  free  space, 
but  with  touching  detail  domain  boundary  as  illustrated  in  Figure  18-19.  Then,  an  antenna  on 
platform  problem  is  simulated  using  fully  overlapping  DDM  and  compared  to  the  HFSS 
simulations  in  Figure  20. 


Excitation 


Figure  18:  The  meshing  of  a  dipole  antenna  with  touching  detail  domain  boundary. 


In  Figure  19(a),  iterative  convergence  is  plotted  as  well  as  the  input  impedance  comparison  with 
HFSS  simulations  at  a  low  frequency.  The  magnitude  of  the  electrical  fields  is  plotted  in  Figure 
19(b).  We  note  that  there  are  some  evanescent  discontinuities  residing  on  the  touching  detail 
domain  boundary,  probably  due  to  the  completely  non-conformal  meshing  in  the  detail  and 
background  domains. 
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(a) 


(b) 


Figure  19:  (a)  iterative  convergence  and  impedance  comparison;  (b)  magnitude  of  the  electrical  fields;  of  the 
dipole  in  free  space,  but  with  touching  detail  domain  boundary. 


In  Figure  20,  an  antenna  placed  on  a  PEC  platform  is  simulated  using  fully  overlapping  DDM 
approaches  at  a  low  frequency  (30  MHz).  The  antenna  input  impedance  and  directivities  in  x—  z 
and  y  —  z  planes  are  compared  to  HFSS  simulations.  Better  accuracy  can  be  expected  by 
increasing  meshing  density  and  radiation  truncation  boundary. 


Corporation  XY  Plot  2  HFSSDesign2 


Zddu  =  2.41  -  226.53/  ZHFSS  =2.65  -  236.7 9j 

Figure  20:  The  directivity  and  input  impedance  comparison  between  fully  overlapping  DDM  and  HFSS 
simulations  in  simulating  a  dipole  on  a  PEC  platform  at  30  MHz. 


In  conclusion,  fully  overlapping  DDM  possesses  the  property  of  faster  convergence,  good 
accuracy  and  supreme  modeling  flexibility.  It  is  naturally  suitable  for  parallel  computing  and 
antenna  in  situ  design  optimizations.  The  accuracy  of  the  approach  is  illustrated  by  the  h- 
refinement  convergence,  which  shows  the  error  of  our  fully  overlapping  DDM  can  be  controlled 
in  almost  exactly  the  same  rate  as  conventional  FEM.  Moreover,  significant  faster  iterative 
convergence  distinguishes  our  approach  from  traditional  non-overlapping  DDM.  Finally,  as  a 
demo  example,  an  antenna  on  a  platform  structure  is  simulated  using  fully  overlapping  DDM  and 
compared  to  the  standard  HFSS  results. 
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Research  Summary 


Modeling  of  Electromagnetic  Properties  of  Magnetoelectric  Thin  Films: 

We  present  some  breakthroughs  in  the  numerical  and  theoretical  modeling  of  magnetoelectric 
media.  Our  initial  investigation  was  the  implementation  of  planar  waveguide  structures  using 
magnetoelectric  media.  This  investigation  was  necessary  to  understand  electromagnetic  wave 
propagation  in  magnetoelectric  media,  and  then  to  help  propose  microwave  devices  based  upon 
the  propagation  phenomenon.  Electromagnetic  wave  propagation  within  magnetoelectric,  or 
bianisotropic  materials  in  general,  have  been  investigated  since  the  late  1960s,  and  though  there 
seems  to  be  a  host  of  potential  microwave  and  millimeter  wave  applications,  there  has  been  little 
progress  in  terms  of  applications  because  of  the  lack  of  progress  in  the  actual  implementation  of 
this  material.  Bianisotropic  media  represents  materials  where  the  electric  and  magnetic  fields  are 
coupled  in  its  constitutive  relations,  as  shown  in  (1)  below. 

D=fE+|H  ^ 

B  =  //H  +  gE 

The  permittivity,  permeability  and  magnetoelectric  parameters  are  usually  3x3  matrices.  This 
fact  makes  it  rather  difficult  to  obtain  analytic  solutions  for  electromagnetic  waves  travelling  in 
planar  structures.  Implementation  of  magnetoelectric  material  by  many  groups,  including  Prof. 
Carmen  at  UCLA,  has  introduced  renewed  interest  in  the  material. 

In  previous  reviews,  we  had  presented  a  numerical  modeling  tool  to  obtain  solutions  to 
electromagnetic  wave  propagation  in  complex  materials.  Using  this  tool,  implemented  using 
COMSOL  multiphysics  software,  we  are  able  to  observe  the  propagation  characteristics  for  a 
magnetoelectric  thin  film  in  a  planar  waveguide  structure.  In  this  research  study,  the  material 
characteristics  of  the  substrate  are  very  important.  We  assumed  the  magnetoelectric  material  will 
have  magnetic  characteristics  as  is  the  case  with  either  the  piezomagnetic  or  magnetostrictive 
materials  used  in  the  realization  of  the  magnetoelectric  effect.  The  permittivity  and  permeability 
of  the  media  used  in  this  study  is  given  as  (DC  magnetic  bias  along  z-axis): 
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Simulation  of  the  planar  bianisotropic  waveguide  is  done  using  lossless  bianisotropic  tensors 
This  implies  that  the  tensors  are  Hermitian  in  form,  and  the  magnetoelectric  tensors  have 
the  relationship  £  =  a  and  g  =a+  The  theoretical  values  for  the  magnetoelectric  tensor  (ME), 


which  represent  the  magnetic  point  group  2,  m’,  or  2/m’,  used  in  this  study  are  given  as: 
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—70.15 
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(3) 


The  above  ME  matrix  has  values  that  are  probably  too  high  and  difficult  to  achieve  with  bulk 
or  layered  materials.  However,  it  is  shown  here  to  investigate  electromagnetic  effects  that  could 
be  achieved  with  these  materials.  We  investigated  a  microstrip  line  as  shown  in  Fig.  la.  This 
planar  structure  consists  of  two  layers  where  the  bottom  layer  is  the  magnetoelectric  medium  and 
the  second  layer  is  a  low  loss  dielectric  (insulator).  The  dielectric  layer  is  used  to  allow  a  better 
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portion  of  the  wave  to  propagate  in  the  dielectric  layer  than  in  the  bianisotropic  layer  where 
attenuation  of  the  wave  may  be  greater.  However,  this  set  up  also  allows  for  the  propagation 
phenomena  existing  in  the  bianisotropic  media  to  show  up  in  the  entire  structure. 

Simulation  of  the  microstrip  line  was  done  at  25GHz.  The  dimensions  of  the  line  are  given  as 
B=0.067X,  D=0.058X,  and  W=0.167L  Figure  lb  shows  the  time  average  power  flow  in  the  z- 
direction  (propagating  out  of  page).  Numerical  analysis  of  the  structure  shows  a  field 
displacement  phenomenon.  This  phenomenon,  although  with  a  DC  magnetic  field  bias  along 
direction  of  propagation,  is  similar  to  that  which  occurs  in  ferrites  and  magnetized  semi¬ 
conductors,  when  the  magnetic  bias  is  transverse  to  the  propagation  direction  and  the  plane  of  the 
thin  films.  Such  field  displacement  effect  is  used  in  several  nonreciprocal  microwave  and 
millimeter  wave  devices. 


Free  space 
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Fig.  1:  Schematic  diagram  of  the  microstrip  line  is  shown  in  (A).  The  red  lines  in  the  figure  represent 
perfect  electric  conductors  (PEC).  (B)  Simulated  result  for  the  magnetoelectric  microstrip  line.  Observe 
propagation  along  left  edge  of  microstrip. 


Fig.  2:  (A)  Time  average  power  flow  across  the  face  of  the  planar  structure.  Maximum  power  density 
along  each  edge  of  the  microstrip  line.  (B)  2D  surface  plot  of  the  propagation  fields  along  one  edge  of  the 
microstrip. 

It  can  be  observed  that  there  is  increased  power  propagation  along  one  edge  of  the  microstrip 
line,  with  much  less  power  propagating  on  the  other  end.  This  magnetoelectric  microstrip  line 
configuration  can  be  easily  applied  to  the  design  of  2-port  isolator  or  a  3-port  circulator.  The 
displaced  power  observed  in  the  structures  stems  from  the  displaced  fields  within  the  structure, 
resulting  from  the  material  properties  of  the  media.  The  magnetoelectric  thin  film  with  a  DC 
magnetic  field  bias  serves  as  a  nonreciprocal  media  causing  the  displacement  of  the  electric  and 
magnetic  field  components  towards  the  right  part  of  the  structure.  A  subdomain  power 
integration  of  each  half  (left  and  right)  of  the  magnetoelectric  microstrip  structure  shows  that 
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99%  of  the  power  is  displaced  towards  the  left  edge  of  the  line.  This  confirms  that  a 
nonreciprocal  field  displacement  effect  is  possible  with  thin  magnetoelectric  films  biased  along 
the  direction  of  propagation. 

An  important  part  of  the  investigation  into  practical  electromagnetic  applications  of  the  ME 
media  is  the  theoretical  model  of  the  media  used  in  numerical  simulations.  Harshe,  [G.  Harshe, 
J.P  Dougherty  and  R.E.  Newnham,  “Theoretical  Modeling  of  Multilayer  Magnetoelectric 
Composites,  Int.  Journal  ofAppl.  Electromag.  in  Mater.  4,  1993 ]  obtained  theoretical  models  for 
the  longitudinal  ME  voltage  coefficient  of  multilayer  composites  composed  of  lead  zirconate 
titanate  (PZT)  and  ferrite  layers  (Cobalt  and  Nickel  Ferrite).  However,  there  were  large 
deviations  between  the  ME  voltage  coefficient  obtained  via  theoretical  modeling,  and 
corresponding  experimental  values  of  fabricated  structures.  The  experimental  results  were 
comparable  to  those  reported  in  literature,  but  the  theoretical  results  were  several  times  higher. 
Reasons  for  the  poor  agreement  between  theoretical  and  experimental  results  include  poor 
interface  coupling  between  the  layers  and  the  application  of  boundary  conditions  used  to  obtain 
the  model.  Since  then,  other  works  have  looked  into  obtaining  theoretical  models  that  account  for 
the  mismatch  at  the  composite’s  bonding  interface.  This  was  done  using  an  interface  coupling 
parameter  that  models  the  strain  transfer  relationship  between  the  piezoelectric  and 
magneto strictive  phases.  The  ME  voltage  coefficients  obtained  in  these  works  had  similar  trends 
as  the  experimental  data,  however,  there  still  remained  large  deviation  in  the  overall  magnitude 
of  the  ME  voltage  coefficients,  as  theoretical  results  remained  several  times  higher  than 
experimental  values. 

In  contrast  to  previous  work,  we  apply  fundamental  electromagnetic  wave  boundary  conditions 
on  the  fields  within  the  composite  structure.  We  introduce  new  theoretical  models  for  the  ME 
effect  in  a  piezoelectric/magnetostrictive  bilayer  that  better  approximates  the  experimental 
results.  The  model  is  obtained  by  solving  the  constitutive  equations  of  each  layer  for  the  all  fields 
present,  and  then  applying  a  field  averaging  method,  along  with  boundary  conditions  on  the 
components  of  the  fields  at  the  composite  interface,  to  obtain  a  homogenized  layer.  The  layer  is 
characterized  in  terms  of  its  effective  permeability,  effective  permittivity  and  the  effective  ME 
susceptibility  tensor  with  constitutive  equations  of  the  form  shown  in  (1). 


We  have  successfully  modeled  the  ME  media  for  each  of  the  three  basic  configurations  as  can  be 
observed  in  Fig.  3.  We  have  improved  the  agreement  with  measured  data  obtained  in  the  work 
done  by  Harshe.  We  summarize  here  the  constitutive  relationship  for  each  of  these  cases.  First, 
for  the  longitudinal  case,  we  have  the  media  poled  and  biased  in  the  3  direction.  The 
homogenized  constitutive  relationship  is  obtained  as 
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(C) 
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Fig.  3:  Orientation  of  biasing  fields  for  the  magnetoelectric  media  in  terms  of  the  poling  (electric)  and 
biasing  (magnetic)  fields.  (A)  Longitudinal;  (B)  In-plane  longitudinal;  (C)  Transverse.  Each  orientation 
produces  a  different  form  of  the  constitutive  relationship. 


Using  the  magnetoelectric  susceptibility  tensor,  we  can  compute  the  magnetoelectric  voltage 
coefficient.  Plots  of  the  calculated  ME  voltage  coefficient  is  shown  in  Fig.  4.  We  also  consider 
the  effects  of  strain/stress  coupling  between  the  piezoelectric  and  magneto strictive  material.  We 
make  use  of  an  interface  coupling  parameter,  k,  to  model  the  strain  transfer  between  phases.  An 
interface  coupling  factor  of  one  (k=l)  implies  perfect  coupling  and  thus  complete  strain  transfer 
between  phases.  It  is  important  to  note  that  one  cannot  expect  a  perfect  coupling  at  the  interface 
since  the  phases  are  not  joined  at  a  microscopic  level.  In  most  cases,  some  form  of  adhesive  is 
used  to  hold  the  phases  together. 


Piezoelectric  Volume  Fraction  Piezoelectric  Volume  Ftadion 


Fig.  4:  Plots  of  the  longitudinal  ME  voltage  coefficients  for  PZT-CFO  bilayer.  (A)  shows  the  effect  of  the 
interface  coupling  parameter  between  phases.  (B)  compares  our  theoretical  model  to  measured  results 
and  previous  theoretical  model  by  Bichurin  et  al. 

The  theoretical  model  here  obtained  better  characterizes  the  ME  voltage  coefficient  of 
piezoelectric  and  magnetostrictive  bilayers.  This  new  theoretical  model  realized  here 
significantly  reduces  the  deviation  between  the  magnitude  of  previous  theoretical  models  and 
experimental  data,  while  maintaining  established  trends  of  the  ME  voltage  coefficient.  The 
measured  value  of  the  ME  media  for  an  equal  volume  composition  of  the  piezoelectric  and 
magnetostrictive  phases  was  74  mV/cmOe  obtained  by  Harshe.  Our  theoretical  model  gives  a 
value  of  140  mV/cmOe  for  perfect  coupling  between  phases  compared  to  340  mV/cmOe  from 
previous  models.  Additional  studies  on  the  coupling  possible  between  phases  can  lead  to  the 
actual  value  of  the  coupling  coefficient,  k,  for  any  two  phases.  From  our  theoretical  model,  we 
can  approximate  the  interface  coupling  between  phases  as  k=0.58.  The  theoretical  model 
includes  expressions  for  the  effective  permeability,  permittivity,  and  ME  susceptibility  tensor  of 
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the  bilayer.  These  expressions  obtained  will  allow  for  better  investigation  of  the  ME  composite 
towards  potential  device  applications.  The  theoretical  model  also  accounts  for  the  bond 
mismatch  at  the  interface  of  the  layers,  which  is  important  in  the  physical  realization  of  an  ME 
composite  with  a  high  ME  voltage  coefficient.  Detailed  information  on  the  theoretical  model  is 
available  in  the  journal  papers  submitted  as  part  of  this  research  effort. 

Theoretical  modeling  for  the  transverse  ME  effect  was  also  completed.  In  this  orientation,  the 
direction  of  the  magnetic  field  bias  is  changed  from  the  3 -direction  to  the  1 -direction,  as 
observed  in  Fig.  3.  The  media  remains  poled  in  the  3-direction.  The  constitutive  relationship  for 
the  transverse  ME  effect  is  given  as 
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Similar  to  the  longitudinal  case,  plots  the  magnetoelectric  voltage  coefficient  for  a  PZT -Nickel 
ferrite  (NFO)  bilayer  can  be  seen  in  Fig.  5. 


Piezoelectric  Volume  Fraction  Piezoelectric  Volume  Fraction 

Fig.  5:  Plots  of  the  transverse  ME  voltage  coefficients  for  a  PZT -NFO  bilayer.  (A)  shows  the  effect  of  the 
interface  coupling  parameter  between  phases.  (B)  compares  our  theoretical  model  to  previous  theoretical 
model  by  Bichurin  et  al.  Note  that  our  model  shows  a  lower  ME  voltage  coefficient  as  is  consistent  with 
measured  values. 


The  results  of  the  modeling  are  consistent  with  the  behavior  of  the  ME  media  from  previous 
theoretical  models.  However,  the  magnitude  of  the  ME  voltage  coefficient  for  our  theoretical 
model  is  much  lower  than  previous  models.  This  is  also  consistent  with  comparisons  of  previous 
theoretical  results  to  measured  data  available  for  the  ME  media.  The  theoretical  model  obtained 
for  the  transverse  ME  effect  includes  expressions  for  the  effective  permeability,  permittivity,  and 
ME  susceptibility  tensor  of  the  bilayer.  From  the  ME  susceptibility  tensor,  the  ME  voltage 
coefficient  is  obtained.  The  transverse  ME  effect  produces  a  higher  ME  susceptibility  compared 
to  that  of  the  longitudinal  case.  This  is  based  solely  on  the  material  property  of  the  media  and 
hence  the  greater  resultant  effect  of  the  magnetic  field  bias  in  the  1  -direction  of  the  media. 

Finally,  we  obtained  a  homogenized  form  for  the  constitutive  relationship  of  the  in-plane 
magnetoelectric  media.  In  this  orientation,  the  media  is  poled  and  biased  in  the  1 -direction  (see 
Fig.  3).  As  was  done  for  each  of  the  ME  orientations,  we  obtained  the  homogenized  effective 
permeability,  permittivity,  and  ME  susceptibility  tensor  of  the  media.  Using  this  information  we 
obtained  the  constitutive  relationship  of  the  media  as  shown  below. 
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In  summary,  we  have  obtained  a  complete  theoretical  model  for  the  magnetoelectric  effect  in 
piezoelectric  and  magneto strictive  bilayers.  The  model  includes  expressions  for  the  homogenized 
effective  permeability,  permittivity,  and  ME  susceptibility  tensor  of  the  bilayer.  We  also  have 
been  able  to  complete  numerical 
modeling  of  a  theoretical 
magnetoelectric  media  using  planar  , 
waveguide  structures.  We  have 
observed  possible  applications 
using  the  non-reciprocal 

propagation  phenomena  of  the 
structure.  However,  the 

magnetoelectric  susceptibility 

obtained  from  the  theoretical  modeling  of  the  PZT-ferrite  bilayers  is  of  small  magnitude  such 
that  similar  effects  in  waveguide  structures  have  not  been  duplicated.  We  are  currently 
investigating  the  use  of  materials  with  higher  piezoelectric  and  piezomagnetic  coupling 
coefficients  to  influence  an  increase  in  the  magnetoelectric  susceptibility  of  the  media. 

Magnetic  Semiconductors  for  Non-reciprocal  Device  Applications: 

Magnetic  semiconductors  have  an  anisotropic  permeability  and  permittivity  under  a  DC  magnetic  field 
bias.  The  shape  of  the  permeability  and  permittivity  tensors  depend  on  the  direction  of  the  DC  magnetic 
field  bias.  Assuming  a  Z-directed  DC  magnetic  field  bias,  the  permeability  and  permittivity  of  the 
magnetic  semiconductor  is  given  by 


Figure  6:  Model  of  magnetic  semiconductor  slotline. 
DC  magnetic  field  bias  is  0.3T  in  the  Y-direction. 
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The  value  of  the  tensor  components  depends  on  the  material  characteristics  of  the  media,  as  well 
the  magnitude  of  the  DC  magnetic  field  bias.  There  are  three  basis  configurations  based  upon  the 
propagation  and  bias  directions:  Faraday,  Perpendicular,  and  Voigt.  Analytic  solutions  are 
obtained  for  all  three  cases  in  unbounded  media.  The  results  obtained  from  the  analytical 
solutions  show  that  magnetic  semiconductors  behave  similar  to  materials  with  solely  anisotropic 
permeability  or  permittivity.  However,  the  propagation  constants  for  the  magnetic  semiconductor 
contain  more  design  variables  that  would  surely  be  useful  as  more  applications  are  discovered. 
The  magnetic  semiconductor  modeled  here  is  a  high  mobility  GaAs/AlGaAs  hetero structure 
doped  with  Manganese.  This  magnetic  semiconductor  was  chosen  for  its  high  mobility,  which 
holds  at  sub-room  temperature  (77K).  Currently,  there  are  not  any  high  mobility  magnetic 
semiconductors  with  Curie  temperature  at  room  temperature.  An  obvious  disadvantage  in  using 
magnetic  semiconductors  or  semiconductors  in  general,  in  planar  waveguide  structures  is  the 
ohmic  dissipation  that  is  introduced  into  the  structure.  For  this  reason,  we  utilize  an  insulating 
low  loss  dielectric  layer  ( sr  =  2.7)  to  help  alleviate  this  problem  (see  Fig  6).  The  dielectric  layer 

is  a  propagation  channel  which  holds  a  majority  of  the  electromagnetic  waves  in  the  structure. 
Using  the  properties  of  the  magnetic  semiconductor  and  the  magnitude  of  the  external  DC 
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magnetic  field  bias  as  0.3T  in  the  Y-direction,  we  obtain  the  permittivity  and  permeability  of  the 
material  as  follows: 
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Simulation  of  the  planar  structure  was  carried  out  at  30GHz  using  COMSOL  multiphysics  solver. 
We  observe  a  displacement  of  the  power  propagating  along  the  slot,  towards  one  edge  (see  Fig. 
7(a)).  The  reason  for  the  observed  field  displacements  comes  from  the  interaction  between  the 
Hall  currents,  produced  by  the  application  of  the  external  DC  magnetic  fields  to  the  magnetic 
semiconductor,  and  the  electric  fields  propagating  in  the  structure.  To  see  if  the  magnetic 
semiconductor  provides  larger  displacement  than  a  non-magnetic  semiconductor,  we  compare 
the  power  displaced  in  the  structure  to  that  of  an  often  used  GaAs  semiconductor.  This 
semiconductor  has  similar  elemental  base  as  the  magnetic  semiconductor  heterostructure.  The 
material  characteristics  of  the  semiconductor  are  given  as:  carrier  concentration  n  =  2x10 20 cm-3 ; 
the  electron  mobility  //,  =  20 in /V  -s  at  77K;  static  dielectric  constant  es0  =13.1;  the  effective 

mass  m*  =  0.063 .  A  comparison  of  the  results  comparing  the  nonreciprocal  phenomenon  in  both 
substrates  is  shown  in  Fig.7.  We  observe  that  there  is  significant  improvement  in  the  displaced 
power. 


(a)  (b) 


Figure  7 :  Comparison  of  the  power  displacement  in  a  magnetic  semiconductor  slotline  to  a  semiconductor 
slotline:  (a)  Magnetic  Semiconductor,  (b)  Semiconductor 

The  phenomena  shown  here  is  readily  applicable  to  non-reciprocal  devices  such  as  isolators  and 
circulators.  A  3D  numerical  model  for  a  3  port  circulator  currently  under  development  is  shown 
Fig.  8. 


Figure  8:  3-port  circulator  design  based  upon  magnetic  semiconductor  slotline,  (a)  3D  simulation  model; 
(b)  excitation  at  port  1;  (c)  excitation  at  port  2;  (d)  excitation  at  port  3 

Characterization  of  EM  Properties  of  On-Wafer  Magnetic/Dielectric  thin-films 
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Error  analysis:  Our  newly  developed  method  for  magnetic  dielectric  on-wafer  thin-film 
characterization  is  discussed  here.  This  method  requires  the  generation  of  various  test  on  the  test 
wafer.  However,  the  generated  fixtures  will  always  have  some  errors  due  to  the  imperfection  of 
the  fabrication  process.  In  our  method,  microstrip  or  other  planar  transmission  lines  can  be  used; 
however,  this  discussion  is  limited  to  microstrip  lines.  Since  the  width  of  the  microstrip  line 
plays  a  very  important  role,  it  is  imperative  to  present  an  error  analysis  in  terms  of  this 
parameter.  Keeping  in  mind  that  two  different  widths  of  the  microstrip  line  are  used  in  our 
method,  we  generated  two  sets  of  random  numbers  for  the  microstrip  line  widths.  In  this 
analysis,  we  used  an  electromagnetic  simulation  tool  to  study  the  errors  in  measuring  the 
permittivity  and  permeability  due  to  errors  in  the  microstrip  line  width.  In  the  simulation,  we 
used  two  microstrip  transmission  lines  which  had  line  widths  of  500  pm  and  600  pm, 
respectively  (see  Fig.  9).  Two  sets  of  random  numbers  were  generated.  The  first  set  had  a  mean 
of  500  pm  and  a  standard  deviation  of  ±2.5  pm.  Similarly,  the  second  set  had  a  mean  of  600  pm 
and  a  standard  deviation  of  ±3  pm.  Fig.  10  shows  simulated  results  with  a  standard  error 
analyses  for  both  s\  and  //  The  red,  blue,  and  green  curves  represent  average  values,  upper 

95%  limit,  and  lower  95%  limit,  respectively.  We  used  standard  error  given  by  SE  =  s/yfn 
where  s  is  the  sample  standard  deviation  and  n  is  the  sample  size.  Also,  the  upper  and  lower  95% 
limits  are  given  by:  95%  confidence  limits  =  x±(S£'x  1.96)  where  x  is  the  sample  mean.  The 

maximum  and  minimum  variations  of  s'r  for  the  upper  and  lower  95%  limits  are  0.182  and  0.180, 
respectively.  Also,  the  maximum  and  minimum  variations  of  /Lt'r  for  the  upper  and  lower  95% 
limits  are  0.205  and  0.2012,  respectively.  The  expected  maximum  and  minimum  errors  of 
extraction  value  of  s'r  with  respect  to  the  exact  value  of  3  in  the  upper  and  lower  95%  confidence 
limits  are  3.73%  and  2.78%,  respectively.  The  expected  extraction  error  with  respect  to  the  exact 
ft',  value  of  2  has  15.13%  for  the  maximum  error  and  4.10%  for  the  minimum  error.  The  results 
of  the  error  analysis  for  the  extraction  of  e'r  and  ju'r  show  that  our  newly  developed  magnetic 
dielectric  characterization  method  can  yield  extracted  values  for  both  e\  and  //,  with 
approximately  10%  error  even  if  the  test  structures  have  fabrication  errors. 


Fig.  9:  Test  set  of  microstrip  lines  with  different  line  widths  (W1  =  500  //m 
and  W2  =  600  fim) 
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Another  consideration  is  a  standard  error  analysis  of  the  material  loss.  Fig.  11  shows  the 
simulated  results  with  standard  error  analysis  for  both  s "r  and  // (the  exact  values  of  e"r  and  n"r 
are  0.015  and  0.01,  respectively).  This  analysis  includes  all  possible  variables  in  both  DUTs  and 
the  TRL  calibration  kits.  Although  this  analysis  includes  all  the  variables,  the  maximum  and 
minimum  variations  for  the  upper  and  lower  95%  confidence  limits  are  very  small.  For  e"r,  the 
maximum  and  minimum  variations  are  0.000426  and  0.000254,  respectively.  Also,  the  maximum 
and  minimum  variations  for  p  "r  are  0.000847  and  0.000495,  respectively.  The  variations  in  // 
are  larger  than  they  were  for  the  s"r  analysis  case;  however,  the  sizes  of  the  variations  were  still 
small. 

Anisotropic  thin-film  characterization:  The  accurate  characterization  of  the  electromagnetic 
properties  (s  and  //)  of  new,  on-wafer  thin  films  is  crucial  for  accessing  their  potential  use  in  the 
design  of  microwave  devices,  antennas,  and  a  variety  of  sensors.  Furthermore,  many  of  the  new 
materials  being  developed  are  anisotropic,  and  conventional  techniques  for  isotropic  material 
characterization  methods  cannot  be  used.  The  method  developed  in  our  research  is  based  on  the 
mapping  of  an  anisotropic  region  in  the  Z-plane  into  an  isotropic  region  in  the  W-plane. 
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Fig.  10:  Simulated  results  for  the  real  parts  of  permittivity  (left)  and  permeability  (right)  with  standard 
error  analysis  (s'r  =3  and  ju'r= 2  are  the  exact  values).  Errors  introduced  in  both  DUT’s  and  TRL 
calibration  kits 
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Fig.  1 1 .  Simulated  results  for  the  imaginary  parts  of  permittivity  (left)  and  permeability  (right) 
with  standard  error  analysis.  Red  curves  on  both  plots  are  the  average  values  of  e"r  and  ju"r.  Blue 
and  Green  curves  represent  the  upper  and  lower  95%  limits,  respectively.  (e"r=0.015  and  ju"r= 0.01 
are  the  exact  values). 


Fig.  12  shows  a  cross  section  of  the  microstrip  in  the  Z-plane  and  the  W-plane.  The  effective 
permittivity  can  be  determined  from  the  measurements  of  microstrip  lines,  however,  the  effective 
permittivity  is  a  function  of  both  sg  and  He  which  are  two  unknowns.  Thus,  it  is  necessary  to 
measure  two  different  microstrip  lines 


Fig.  12:  Cross  section  of  (a)  microstrip  on  123456739  10 

isotropic  substrate  Frequency  (GHz) 


Fig.  14.  Simulation  result  for  biaxial  anisotropic  substrate 
to  find  two  different  effective  (£xt=3,  %=6  and  %=9) 

permittivity  values,  then  it  is  possible  to 

solve  the  unknowns,  such  as  sg  and  He.  From  the  determined  eg  and  He,  the  anisotropic  material 
permittivity  can  be  found. 
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sx  =  s  ^  £_  for  uniaxial  anisotropic 
sx  ^  sy  ^  s_  for  biaxial  anisotropic 


(9) 


Since  the  optical  axes  of  the  anisotropic  material  (not  known)  and  the  measurement  axes  of  the 
measurement  system  are  not  the  same,  the  matrix  will  be  in  general  non-diagonal.  Fig.  13  shows 
the  angles  between  the  xyz  and  the  x'y'z'  system  where  9  is  the  rotation  angle  along  the  z-axis 
and  that  (p  is  the  rotation  angle  of  the  x-axis. 

For  uniaxial  anisotropic  material,  the  permittivity  tensor  can  be  found  from  the  measured 
effective  permittivity  of  two  different  microstrip  lines  because  in-plane  permittivity  tensor 
elements  are  same  (ex=£y).  However,  it  is  necessary  to  consider  two  different  propagation 
directions  on  the  in-plane  (x-y  plane)  optical  axes.  One  is  the  x-axis  propagation,  and  another  is 
the  y-axis  propagation.  Each  propagation  direction  can  be  considered  as  a  microstrip  line  on  a 
uniaxial  anisotropic  substrate  problem. 
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Fig.  13.  Misalignment  angle^g  jg  Simulated  results  for  an  anisotropic  substrate^ z 
system)  and  the  measuremeincluding  ^  (ex_3j  %=6,  ^-9,  and  ^=1.5  are  the  actual,sor 
(diagonal  elements  only)  in  thValues  in  the  simulation  :ero 


We  tested  our  method  using  electromagnetic  simulation  tools  with  arbitrary  biaxial  anisotropic 
substrate  with  sx= 3,  sy= 3,  and  ez=9.  In  this  simulation,  it  is  assumed  that  the  optical  axes  of  the 
biaxial  anisotropic  substrate  are  known.  Fig.  14  shows  the  simulated  result  for  characterization  of 
biaxial  anisotropic  substrate.  The  variation  of  extracted  sx  was  from  3.044  to  3.081  over  the 
frequency  range  of  1  GHz  to  10  GHz.  Also,  the  extracted  values  of  sy  and  s-  varied  from  5.770  to 
5.832  and  9.059  to  9.105,  respectively.  All  the  extraction  results  were  below  5%  of  the  relative 
error. 


We  discussed  both  uniaxial  and  biaxial  anisotropic  material  characterization  using  microstrip 
lines.  The  main  assumption  of  the  previous  analysis  was  that  we  knew  the  optical  axes  of 
anisotropic  materials.  If  we  still  assume  that  the  z-axis  is  one  optical  axis,  the  permittivity  tensor 
will  have  non-zero  off-diagonal  elements  and  the  permittivity  tensor  is  given  by: 
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We  simulated  test  structures  on  an  anisotropic  substrate  that  had  sxx  =  3,  syy  =  6,  s--  =  9,  and  sxy  = 
1.5  with  a  thickness  of  100  pm.  Fig.  15  shows  the  extracted  results  for  sxx,  syy,  szz,  and  sxy.  The 
extracted  value  of  sxx  varied  from  3.025  to  3.047  over  the  frequency  range  of  1  GHz  to  10  GHz. 
The  minimum  and  maximum  errors  were  0.846%  and  1.571%,  respectively.  The  extracted  result 
for  syy  varied  from  5.817  to  5.792,  and  szz  varied  from  9.056  to  9.094.  Let  us  consider  the 
simulated  result  for  off-diagonal  element,  sxy.  The  extracted  value  of  sxy  varied  from  1.406  to 
1.407  over  the  frequency  range  of  1  GHz  to  10  GHz.  The  minimum  and  maximum  errors  were 
6.229%  and  6.266%,  respectively. 


On-wafer  characterization  of  sapphire  film:  A  good  example  of  anisotropic  thin-film  material 
is  a  sapphire  wafer.  The  dielectric  constants  of  sapphire  are  1 1.6  for  the  parallel  to  the  c-axis  and 
9.4  for  the  perpendicular  to  the  c-axis.  Fig.  16  shows  a  conventional  unit  cell  of  a  single  crystal 
sapphire  with  the  orientation  of  C-plane  and  i?-plane.  The  permittivity  tensor  of  C-plane  sapphire 
wafer  is  given  by: 
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exx,  syy,  and  ezz 


(11) 
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The  permittivity  tensor  of  C-plane  sapphire  wafer  has  same  form  of  uniaxial  anisotropic 
permittivity  tensor.  The  permittivity  tensor  of  /?-plane  sapphire  wafer  can  be  easily  calculated  to 
be: 
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(12) 


Although,  we  know  the  permittivity  tensor  of  /?-plane  sapphire  theoretically,  it  is  impossible  to 
build  test  structures  on  the  wafer  perfectly  aligned  with  the  optical  axes  and  the  measurement 
axes.  Thus,  there  is  a  misalignment  between  the  optical  axes  and  the  measurement  axes.  As  a 
result,  the  permittivity  tensor  of  /?-plan  sapphire  wafer  will  have  non- zero  off-diagonal  elements. 
For  example,  if  a  microstrip  line  lies  on  /?-plane  sapphire  wafer  with  in-plane  misalignment 
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Fig.  17.  Simulated  results  for  /?-plane  sapphire  wafer  with  misalignment  angles:  sxx= 9.6801, 
£^=10.6882,  £zz=10.0316,  sxy=- 0.6007,  sxz=0A206,  and  £^=-0.9021  are  the  actual  values 


angle  of  25°,  the  permittivity  tensor  becomes: 
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We  tested  our  method  for  characterization  of 
anisotropic  material  (sapphire)  with  misalignment 
angles  using  electromagnetic  simulation  tools.  In 
this  simulation,  we  used  the  permittivity  tensor  in 
(13)  and  Fig.  17  shows  the  extracted  results  for  all 
the  tensor  elements.  The  extracted  value  of  sxx 
varied  from  8.977  to  8.817  over  the  frequency 
range  of  1  GHz  to  10  GHz.  The  extracted  result  for 
£yy  varied  from  9.941  to  9.976  and  szz  varied  from 
10.167  to  10.252.  The  relative  errors  for  the 
extracted  results  of  diagonal  elements  are  lower 


Fig.  16.  Orientation  of  C-plane  and  R-plane  in 
the  conventional  unit  cell  of  a  single  crystal 
sapphire  ( a ,  b,  and  c  are  the  optical  axes  of 
sapphire  crystal) 
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than  9%.  The  extracted  results  for  off-diagonal  elements  were  also  close  to  the  exact  values  in 
the  simulation. 
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Research  Summary 

The  following  major  research  breakthroughs  were  made  in  this  project. 


1.  The  method  of  the  production  of  conducting  textiles  was  developed  (Figure  1). 

2.  Conductivity  of  the  SWNT  composites  reached  5*105  S/cm. 

3.  High  frequency  conductivity  was  evaluated  and  found  identical  to  that  at  low  frequencies. 

4.  Direct  write  approach  to  antenna  preparation  was  achieved  (Fig.2) 

5.  Textiles  and  fibers  with  the  conductivity  of  copper  were  developed. 
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We  had  3  major  accomplishments: 


1) We  quantitatively  mapped  the 
theoretical  limits  of  multi¬ 
functional  performance  (strength- 
to-weight  versus  electrical 
conductivity)  of  the  e-fibers 
relevant  for  UAV  structural 
antenna  applications 

2)  We  produced  copper-coated,  multi¬ 
strand  kevlar  fibers  exhibiting 
extremely  high  strength-to-weight 
ratios  in  combination  with  metallic 
conductivity 

3)  We  successfully  realized  for  the 
first  time  a  flexible,  stretchable, 
woven  broadband  antenna  using 
these  e-fibers. 
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Our  calculations  (Fig.  MSI)  show  that  such  multi¬ 
filament  geometries  should  be  vastly  superior  to 
conventional  all-metal  wires  when  multi¬ 
functionality  is  desired,  such  as  in  structurally 
integrated  antennas  for  small  UAV  applications. 
Copper  coated  on  multi-strand  Kevlar  fibers  yields 
excellent  performance  (i.e.  sufficiently  high 
conductivity  for  an  antenna,  and  retains  the  high 
strength-to-weight  ratio  of  the  substrate  fiber), 
while  use  of  higher  strength  substrates,  such  as 
Zylon  is  expected  to  yield  even  better  performance. 
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Using  copper-coated  kevlar  fibers,  we  demonstrated  for 
the  first  time  a  woven,  flexible  and  stretchable, 
broadband  antenna.  We  quantified  the  effect  of 
weaving  on  the  electrical  and  mechanical 
characteristics  of  the  e-fibers  (Fig.  MS2).  We 
discovered  that  considerable  mechanical  stress  placed 
on  the  e-fiber  during  the  weaving  process  can  result  in  f 
a  decrease  of  electrical  conductivity,  as  well  as 
breakage.  We  developed  an  encapsulating  strategy  to 
maintain  fiber  integrity  throughout  the  weaving 
process,  and  also  showed  that  the  fiber  can  be  repaired 
in  situ  if  breakages  occur. 

Subsequently,  when  the  woven  antennas  were  tested  by  OSU,  the  frequency  response  (Fig.  MS3) 
of  the  broadband  antenna  clearly  showed  the  great  potential  of  the  multi-material  fiber 
architecture,  expected  to  yield  vastly  improved  performance  over  conventional  conductors  for 
load-bearing  applications. 

This  work  resulted  in  a  poster  presentation  at  the  Fall  2009  MRS  meeting  (poster  was  nominated 
for  an  award),  and  a  joint  publication  with  the  Bayram  and  Volakis  groups  at  OSU  (in 
preparation). 
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RESEARCH  SUMMARY 

In  this  project,  we  have  developed  a  coupled  electromagnetic-mechanical  finite  element 
model  that  incorporates  a  large  deformation  low  frequency  mechanical  vibration  module  in  the 
ultra-high  frequency  (UHF)  electromagnetic  field  response.  In  this  simulation,  a  staggered 
coupling  scheme  is  applied.  In  each  time  step,  the  mechanical  displacement,  which  is  assumed  to 
be  independent  to  electromagnetic  field  by  ignoring  the  Lorentz  force,  is  solved  first.  Then  the 
electromagnetic  field,  which  is  influenced  by  mechanical  deformation,  is  solved  in  total 
Lagrangian  method  to  overcome  complicated  boundary  shape.  This  scheme  can  obtain  accurate 
results  and  also  eliminate  unnecessary  sacrifice  in  calculation.  In  this  case,  the  electromagnetic 
field  in  free  space  in  which  the  hyper-elastic  conductor  is  embedded  is  very  important.  But  the 
mesh  in  free  space  can  be  entanglement  due  to  the  large  deformation.  So  a  mesh  adaptive 
method,  R-method  is  proposed  in  this  project. 

The  coupling  of  the  two  disparate  fields  (mechanical  and  electro -magnetic)  has  a  physical 
limitation  for  UHF,  since  the  electro -magnetic  wavelength  from  UHF  is  much  shorter  than  the 
antenna  constitution.  As  such,  the  structural  dimension  cannot  be  traversed  in  a  period  of 
electromagnetic  field.  An  important  objective  that  has  been  executed  throughout  the  research  is 
to  propose  a  multi-time  scale  modeling  for  physical  phenomena  at  two  very  different  driving 
frequencies  that  result  in  a  significant  scale  difference.  The  multi-time  scale  models  being 
developed  is  proposed  to  resolve  this  frequency  difference.  Various  modules  developed  in  this 
task  are  outlined  below. 

1.  Computational  Modeling  of  the  Electro-Magnetic  Fields  using  a  Weak  Form  of  the 
Maxwell’s  equation 

The  computational  modeling  contains  two  dependent  fields,  mechanical  vibration  and 
electromagnetic  response.  To  better  illustrate  the  idea  of  load  bearing  antenna,  large  deformation 
assumption  is  considered  in  the  mechanical  field  part.  We  apply  Saint  Venant-Kirchhoff  model 
as  the  hyperelastic  constitutive  model 
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Here  S_  is  2nd  Piola-Kirchhoff  stress  tensor,  E  is  Lagrangian  Green  strain  tensor,  W  is  strain 
energy  density  function,  u  is  the  displacement,  A.  and  \i  are  Lame  constraints. 

Construct  potential  energy  of  the  mechanical  field  of  the  system  as. 


]!(«)=  |  w(c(u))dV-  J  u  pobdV-  J  ut_0dS 

n0  dn„ 


(1.5) 


Here  t0  is  the  traction,  b  is  the  body  force. 

Then  do  the  variation  of  it,  obtain  the  weak  form  of  the  mechanical  field, 

^n(w)  =  J  (Su  p0v)dV+\[slL  -pyv-^du  pab)fN-  J  5uTodS=0  (1.6) 

n„  f>„  f>„  fii0 

With  the  weak  form,  we  discretize  the  domain.  Fixed  the  right  edge  and  apply  a  sinusoidal 
loading,  as  shown  in  Figure  1 . 


Figure  1.  Mechanical  excitation  on  conductor’s  edge 

We  note  that  convergence  is  related  to  the  number  of  elements.  Figure  2  gives  the  tip 
displacement  convergence  information.  In  this  simulation,  the  most  important  part  is  to  assign 
time  dependent  electromagnetic  field  according  to  mechanical  excitation.  To  accomplish 
commercial  software’s  incapability  in  solving  multi-physics  simulation,  we  have  developed  the 
coupled  electromagnetic  and  mechanical  finite  element  model  not  only  with  quasi-static 
assumption  but  also  the  transient  problem. 
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Tip  Displcement 

Figure  2.  Tip  displacement  infinite  element  model  with  different  number  of  elements 
We  start  with  Maxwell  equation, 


V-(6(x,t))  =  0 

(1.7) 

V-(d(x,t))  =  qe 

(1.8) 

Vxh(x,t)  =  d-^  +  j(x,t) 

(1.9) 

Vxe(x,,)  =  -8bM.') 

-v-  ;  dt 

(1.10) 

where  d  is  electric  flux  density  vector,  e  is  electric  field  intensity  vector,  b  is  magnetic  flux 
density  vector,  h  is  magnetic  field  intensity  vector,  j  is  electric  current  density  vector  and  qe  is 
electric  charge  density.  Since  the  model  is  described  as  a  load  bearing  antenna,  the  boundary  of 
the  structure  would  change  significantly  during  deformation.  Thus,  to  achieve  an  efficient 
algorithm,  total  Lagrangian  method  is  adopted.  In  this  case,  the  configuration  employed  to 
calculate  the  electromagnetic  field  is  always  reference  configuration.  In  order  to  connect  current 
configuration,  which  is  the  deformed  configuration,  to  reference  configuration  is  Nanson’s 
formula  as  shown  following, 

dv  =  JdV 

(1.11) 

nds  =  j(N*F~1)dS 

(1.12) 
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J  =  det(F) 


(1.13) 


Here  the  small  characters  stand  for  the  current  configuration,  and  the  capital  characters  stand 
for  the  reference  configuration,  n  is  the  normal  vector,  v  and  5  are  the  volume  and  surface  of  the 
domain,  separately.  Integrate  the  Maxwell’s  equations  and  apply  the  Nanson’s  formula,  we 
obtain  a  set  of  Maxwell’s  equations  of  reference  configuration, 

V(F(X,t))  =  0  (1.14) 

V(D{X,t))  =  Qe  (1.15) 

VxH(X,t)  =  +  (1.16) 

Vx£(l,t)  =  -® l!l  (1.17) 

Capital  characters  indicate  the  variables  in  reference  configuration.  Although  there  are  no  formal 
definitions,  the  relationship  between  the  variables  of  current  configuration  and  reference 
configuration  are. 


B(X,t)  =  JF~lb(x,t) 


D(X,t)  =  JF~ld(x,t ) 
E(X,t )  =  [  e(x(X,t),t)  +  d-^^ 


(1.18) 

(1.19) 


F  (1.20) 


The  electromagnetic  field  responses  affected  by  mechanical  vibrations  is  observed  by  this 
simulation.  For  UHF  case,  the  wavelength  of  the  electromagnetic  field  is  much  shorter  than 
antenna  structural  dimension.  One  period  of  electromagnetic  field  of  UHF  does  not  cover  a 
whole  structure.  UHF  electromagnetic  field  varies  within  an  antenna  structure.  This  means  that 
the  quasi-static  assumption  which  is  usually  applied  to  electric  mechanical  device  should  not  be 
considered  for  this  particular  application.  Our  simulation  assumption  is  that  (1)  mechanical 
vibration  and  electromagnetic  field  are  not  history  dependent,  and  (2)  conductor  material  does 
not  have  memory,  that  is  non-dielectric,  and  has  no  surface  current  and  charge  density. 
Moreover,  to  simplify  the  electromagnetic  response,  vector  potential  A  and  scalar  potential  cp  are 
introduced.  They  are  defined  as. 


E  = 


(1.18) 
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B  =  VxA 


(1.19) 


If  we  also  assume  that  no  surface  charge  density  exists,  Maxwell’s  equations  in  reference 
configuration  become. 
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H(X,t)  = 


( V  {-VO  -  A  -  (F_1  •  d)  x  (V  x  A)}  •  C"1  Wf_1  •  w)  + -!- {(V  x  A)  •  C} 
1  Fo^ 


(1.21) 


With  these  two  governing  equations  of  electromagnetic  field  coupled  with  mechanical  excitation, 
the  weak  form  si  give  by 
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Here  £0  is  free  space  permittivity,  ju0  is  free  space  permeability,  F  is  the  deformation 
gradient,  C  is  the  right  Cauchy  strain  tensor,  and  j  is  det(F) ,  p0  is  material  density,  b0  is  a 
gravitational  body  force,  /  is  Lorentz  force,  f  is  a  surface  traction.  And  dO  and  S A  are  arbitrary 
fields. 

Along  with  the  variation  with  respect  to  displacement  as  equation  (1.6),  the  coupled 
mechanism  is  described  completely. 

To  validate  the  governing  equation  of  electromagnetic  field,  one  can  appeal  to  commercial 
software  by  switch  off  the  mechanical  excitation.  Although  it  cannot  support  transient  analysis, 
the  stable  value  is  also  useful  to  compare.  Figure  3  shows  the  comparison  of  our  codes’  results 
and  COMSOL’s  result.  The  variable  we  are  interested  in  are  the  magnetic  and  electric  field 
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X  direction  (A/m) 


generated  by  constant  current  input  in  certain  direction.  Note  the  DOF  also  influence  the 
convergences  as  we  discussed  in  mechanical  field. 


Time  (sec)  Time  (sec) 


Figure  3.  Magnetic  and  electric  field  generated  by  constant  current  under  sinusoidal  vibration  results 
and  comparison  with  COMSOL 

From  Figure  3,  we  note  that  the  result  is  much  closer  to  that  of  the  commercial  codes.  And  we  can 
add  the  response  from  mechanical  excitation  which  COMSOL  cannot  fully  compatible.  Till  now, 
the  problem  of  coupled  electromagnetic  and  mechanical  field  simulation  is  complete  in  algorithm. 


2.  Modeling  Free  Space  Domain  Using  R-method  for  Node  Relocation 

The  load  bearing  antenna  on  an  air  foil  or  body  of  aerial  undergoes  mechanical  vibrations 
caused  by  external  excitation  which  moves  in  free  space,  as  shown  in  Figure  4, 


Figure  4.  Antenna  structure  swims  in  free  space 
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As  we  noted  before,  the  governing  equation  for  mechanical  vibration  is  in  the  Lagrangian  reference 
system,  while  the  Maxwell  equations  for  the  electro-magnetic  held  are  in  Eulerian  reference 
system.  Directly  coupling  these  different  frameworks  for  coupled  transient  responses  creates 
inappropriate  computational  domains.  If  governing  equations  in  the  Lagrangian  and  Eulerian 
reference  systems  are  coupled,  inaccurate  data  transfer  can  occur  across  the  surfaces  of  the 
conducting  material  and  free  space.  Thus,  we  have  transformed  Maxwell’s  equations  in  the  format 
of  Eulerian  to  Lagrangian  type.  This  causes  mesh  entanglement  when  the  mesh  of  the  deforming 
antenna  structure  interferes  with  the  static  grid  for  the  free  space.  This  mesh  distortion  and 
entanglement  leads  to  a  negative  Jacobian  in  the  finite  element  transformation.  This  causes  a 
failure  in  computations.  To  overcome  this  problem,  a  local  node  relocation  method  or  r-method, 
is  applied  to  retain  mesh  integrity.  After  one  time  step  of  computation  for  mechanical  deformation, 
nodes  around  surfaces  of  conducting  material  are  re-located  according  to  how  much  nodes  on 
surfaces  of  conductor  move.  Then,  for  newly  located  nodes  of  free  space,  nodal  interpolation  is 
conducted  for  scalar  and  vector  potential  variables  for  electromagnetic  field  computation.  This 
node  relocation  scheme  is  shown  in  Figure  5. 


Figure  5.  Nodal  re-location  and  interpolation:  Newly  located  node  (Solid  circle)  is  interpolated  by 
previous  nodal  information. 


3.  Multi-Time  Scale  Modeling  for  Coupling  Mechanical  and  Electro-Magnetic  Fields 

An  important  part  of  this  research  is  multi-time  scale  modeling  for  coupling  different  physical 
phenomena  governed  by  disparate  frequency  responses.  It  is  prohibitive  to  simulate  mechanical 
and  electromagnetic  time-responses  in  ultra  high  frequency  scale  for  a  finite  element  analysis 
modeling.  To  overcome  this  computational  difficulty,  the  multi-time  scale  modeling  for 
electromagnetic  fields  is  proposed  to  separate  single  time  response  into  time -averaged  response,  a 
function  of  coarse  time  scale,  and  fine  time  response,  a  function  of  coarse  time  and  fine  time 
scales.  Time-averaged  coarse  time  scale  and  fine  time  scale  responses  for  electromagnetic  field  are 
simulated  in  each  time  scale  respectively  with  a  staggered  method. 

The  idea  of  the  multi-time  scale  modeling  is  to  separate  a  variable,  a  function  of  t  (coarse  time) 
and  t  (fine  time),  into  two  variables.  One  is  a  function  of  only  t  and  may  be  periodic  in  a  coarse 
time  scale.  The  other  is  a  function  of  t  and  t  and  almost  periodic  in  the  fine  time  scale.  Applying 
variable  separations  and  respective  assumptions  on  a  governing  equation  separates  the  governing 
equation  into  two  or  more  in  accordance  with  an  increased  number  of  variables.  After  time 
averaging  over  coarse  time  scale,  the  governing  equation  becomes  a  coarse  time  response  equation 
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driven  by  averaged  fine  time  responses.  As  an  example,  suppose  that  governing  equation  is  the 
second  time  derivative  equation. 

ax (t)-y  +  a2 (t)-  y  +  a3(t)-  y  =  f  (; t )  (3.1) 

Define  a  scaling  factor,  g ,  a  variable,  y,  time  derivatives  and  a  time  averaging  operator  <  >,  the 
relationship  between  them  are, 


t  =  g-r,gU  1,  y£ (t,r)  =  M(t)  +  y(t,r)  (3.2) 


The  corresponding  inter-scale  time  derivatives  are  given  as, 

(3.3) 


dy  \  dy 
y  =  —  + — — 
dt  g  dr 


..  d2y  2  d2y  1  d2y 
y  =  — br  + - — +  - 


(3.4) 


dt2  g  dtdr  g2  dr2 
( . ys{t ,  rfj  =  (M  (0)  +  (  y  (t,  t)}  =  M(t) + (  y  (t,  r)j  (3.5) 


Now,  applying  equations  (5)  -  (8)  to  (4),  (4)  yields  an  averaged  governing  equation  as, 


,  .  d2M(t)  ,  .  dM(t)  ......  .  .  d{y{(’T))  .  .  /  ,  o 

ax  (t)  ■  —2  +  a2  ( t )  — +  a3  ( t )  • M(t)  =  -a2  (t)  ■  - — L  -a3(t)-{y  (t,  t)J 


dr 


dt 


dt 


(3.6) 


Equation  (3.6)  is  a  coarse  time  response  equation  and  is  simulated  with  a  fine  time  response  in 
staggered  method.  We  are  currently  working  on  the  implementation  of  this  algorithm  for  the 
coupled  EM-ME  problem. 


4.  Code  Parallelization  for  High  Performance  Computing 

To  carry  out  the  FEM,  an  open  source  academic  code,  FEAP  is  utilized  here.  Based  on  this 
platform,  we  developed  our  own  code  to  achieve  this  project.  Since  the  coupled  electromagnetic 
and  mechanical  finite  element  modeling  contains  conducting  material  domain  as  well  as  free  space 
domain,  it  is  essential  to  construct  a  code  parallelization  to  speed  up  a  simulation  running  time. 
For  mechanical  finite  element  model,  a  large  deformation  model  is  considered  with  Periodic 
vibrations  and  non-iterative  solver  is  used  to  solve  the  transient  problem.  As  the  output  from  the 
Mechanical  Domain  is  the  input  to  the  Electro-Magnetic  Domain,  it  is  extremely  essential  for  us  to 
adopt  domain  decomposition  techniques  in  Mechanical  part  too.  The  following  paragraphs  deals 
with  the  procedures  adopted  to  implement  the  parallel  code  in  Mechanical  Domain  for  a  better 
performance. 

To  develop  the  parallel  implement,  a  serial  input  file  initially  provided  to  FEAP  is  sub-divided 
into  desired  number  of  partitions  depending  on  the  number  of  processors  to  be  used.  For 
example,  if  we  wish  to  run  the  code  on  4  processors  then  we  need  to  divide  the  entire  domain  into 
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four  parts  so  that  each  processor  will  have  its  separate  input  hie.  (Par) METIS  has  been  used  for 
efficient  domain  decomposition  as  shown  in  schematic  diagram. 

Then  we  construct  the  solution  algorithm  for  wide  variety  of  applications.  PETSC,  a 
mathematical  algorithm  aided  by  MPI  for  communicating  information  between  different 
processors  for  efficient  performance  is  appealed.  We  have  also  used  Prometheus  which  can  be 
integrated  easily  with  PETSC  libraries,  as  it  resulted  in  extensive  speed  up.  At  last,  the  graphical 
and  numerical  outputs  of  the  results  are  executed  in  parallel,  either.  Figure  6  shows  the  schematic 
representation  of  the  parallel  computation  with  FEAP. 
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Figure  6.  Schematic  representation  of  the  parallel  computation  with  FEAP. 

After  implementing  the  parallel  code,  we  could  successfully  analyze  the  problem  for  10000 
elements  where  as  the  earlier  serial  version  used  only  225  elements  and  was  comparatively  slower. 
The  results  clearly  showed  improvement  with  the  increase  in  number  of  processors.  The  result 
below  displays  the  time  taken  by  the  processors  in  solving  the  Mechanical  Domain  for  the  first 
1000  time  steps. 


Performance  curve  for  10000  elements 


Figure  7.  Performance  of  certain  model  taken  out  by  different  number  of processors 

We  were  also  successful  in  coupling  the  Mechanical  and  Electro-Magnetic  held  on  a  common 
platform.  The  performance  curve  for  coupled  problem  for  10000  elements  in  the  first  10  steps  is 
shown  below.  We  are  currently  working  on  the  Electro-Magnetic  code  to  improve  its  performance. 
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PERFORMANCE  CURVE  FOR  EM-FEM  CODE 


Figure  8 .  Performance  of EM-FEM  model  taken  out  by  different  number  of processors 


For  validation  purpose,  we  had  selected  a  particular  node  number  and  tried  to  compare  the 
displacements  in  the  Z-direction  for  the  first  five  time  steps.  The  results  proved  to  match  exactly  as 
shown  below. 
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Figure  9.  Finite  Element  simulation  taken  out  by  different  number  of  processors 
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Research  Summary! 

Load-bearing  antennas  are  multi-functional  sensing  (actuating)  and  receiving  (detecting) 
devices  that  are  integrated  with  a  load-bearing  structure.  These  antennas  are  highly 
appealing  for  military  applications,  importantly  Unmanned  Aerial  Vehicles.  In  the 
present  work,  we  have  developed  analytical  techniques  and  computational  tools  for  multi¬ 
scale,  multi-physics  modeling  of  composite  load-bearing  antennas.  The  antenna  structure 
is  subjected  to  mechanical  forces,  temperature  gradients,  and  electromagnetic  fields, 
giving  rise  to  highly  coupled  non-linear  electro-magneto-thermo-mechanical  (EMTM) 
behavior.  Our  model,  based  predominantly  on  first  principles,  employs  the  thermo¬ 
mechanical  governing  equations  (i.e.,  conservation  of  mass,  momentum,  angular 
momentum,  energy  and  second  law  of  thermodynamics)  coupled  with  Maxwell’s 
equations.  The  macro-scale  coupling  terms  for  the  thermo-mechanical  and  electro¬ 
magnetic  interactions  are  deduced  from  micro-scale  behavior.  Our  modeling  has 
identified  92  non-dimensional  numbers  which  quantify  the  competition  between  physical 
effects  in  the  operation  of  load-bearing  antenna.  Depending  on  the  design  of  the  structure 
and  nature  of  the  excitation,  only  a  subset  of  physical  effects  will  be  dominant,  which 
dictates  the  appropriate  computational  model.  A  fixed  relative  ordering  among  all 
competing  effects  determines  a  “regime”  of  antenna/environment  interaction.  The 
mathematical  structure  of  leading  order  equations  for  various  physical  regimes  is 


presented  for  use  in  analysis  and  optimization  of  the  coupled  electro-magneto-thermo- 
mechanical  behavior  of  load  bearing  antennas.  Catalog  of  practical  regimes  of  behavior 
and  corresponding  model  equations  are  presented  and  results  are  analyzed. 


Research  Accomplishments: 

Our  model  employs  the  thermo-mechanical  governing  equations  (i.e.,  conservation  of 
mass,  linear  momentum,  angular  momentum,  energy,  and  the  second  law  of 
thermodynamics)  coupled  with  Maxwell’s  equations.  The  macro-scale  coupling  terms 
for  the  thermo -mechanical  and  electro-magnetic  interactions  are  deduced  from  micro¬ 
scale  behavior.  Our  first-principle  modeling  has  identified  92  non-dimensional  numbers 
which  quantify  the  competition  between  competing  physical  effects  in  the  operation  of 
the  load-bearing  antenna.  Based  on  antenna  structure  and  environment,  we  have  deduced 
a  framework  to  optimize  the  coupled  behavior  and  utilized  this  framework  to  present  a 
catalog  of  all  practical  regimes.  The  mathematical  structures  of  most  relevant  regimes  are 
deduced  and  solutions  are  analyzed. 

Governing  equations  have  been  presented  in  a  non-dimensional  form,  for  example 


Ampere’s  law: 
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The  characteristic  scales  corresponding  to  the  physical  quantity  are  defined  in  table.3  of 
[6],  For  example, 


Charecteristic  Scales 

Notation 

Time  scale  of  electromagnetic  phenomena 

teQ 

Time  scaJe  of  thermo-mechanical  phenomena 

Length  scale  iti  x-direction 

All  the  non-dimensional  qunatities  involved  in  the  coupled  problem  are  listed  below: 
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Table  1:  Non-Dimensional  Quantities 

We  now  define  a  process  and  design  for  modeling  of  load-bearing  antenna  . 


Figure  1:  microsnip  parch  a-ureima 


Designs  of  load-bearing  antenna  are 
largely  microstrip  antennas.  We  model 
a  single  layer  patch  antenna  and  it  is 
further  extended  to  composite  layers 
using  hybrid  methods. 


The  following  restrictions  on  coupled  electro-magneto-thermo-mechanical  behavior 
defines  the  process: 

a)  A  description  of  microstrip  patch  antenna  is  shown  in  Fig.l.  The  antenna  consists 
of  a  dielectric  material  bounded  by  a  ground  plane  at  one  end  and  a  metal  surface 
patch  antenna  with  feedline  at  the  other  end.  The  patch  antenna  surface  is  on  x-y 
plane  (on  the  wing  of  airplane).  We  assume  that  the  antenna  dimensions  along  x 
and  y  »  dimensions  along  z. 

b)  Major  duration  of  flight  of  an  air-plane  is  the  cruise  mode,  where  velocity  of  the 
air-plane  is  constant.We  assume  that  this  velocity  is  directed  along  y-axis  and  of 
the  order  of  lOOm/s. 

c)  Since  velocity  of  air-plane  is  dominant  along  the  direction  of  the  flight,  other 
components  of  velocity  caused  by  deformation  are  of  much  smaller  scales. 
Aerodynamic  thrust  and  weight  of  the  air-plane  are  of  the  order  of  10000  N.  For  a 
small  UAV  the  scales  are  smaller. 

d)  In  defining  scales  for  the  electro -magnetic  problem  it  is  important  to  note  that  the 
fields  vary  drastically  over  the  entire  domain,  so  in  regime  1,  we  consider  only  the 
fields  inside  the  antenna  structure  as  they  contribute  to  the  mechanical  forces.  Far- 
fields  created  by  the  antenna  and  the  distortions  caused  by  antenna  structure  are 
studied  in  a  different  regime. 


56 


e)  The  multi-scale  coupled  behavior  of  an  LBA  requires  that  the  time  scales  are 
different  for  electromagnetic  and  mechanical  problems,  since  EM  fields  vary  at 
much  faster  rates  than  mechanical  fields. 

Typical  values  for  charecteristic  scales  corresponding  to  this  process  are  deduced  as 
shown  in  table  2: 


Table  2:  Non-Dimensional  Quantities 


These  scales  provide  a  relative  ordering  of  the  non-dimensional  quantities  listed  in  table. 

£  =  i±  =  0.01 


1  in  terms  of  slenderness  parameter  ( 
process  is  defined  as: 


)  and  the  regime  corresponding  to  this 


A t  j4ie°  ,  £?3i  =  ,  Bj2  =  Bat 11  ,  Di  =  Die3  ,  jDj  =  D2E3 

F-2,  =  F2Je-*  ,  Fl2  =  FllE-1  ,  =  UZE 3  ,  J3  =  J3e~2 


Su  =  Si  if'1  -.Sn  =  =  S12E0  ,Sl3  =  S3 1  =  SuE3  .S33  =  f>33f2 

Ql  =  Ql  =  QiE  1  1  Pi  =  Pa  =  PjE  5  <  Pij  =  F-ijE  J 


where,  e  <  AL..  Bsl,  £?32,  Du  Da,  F?i ,  F12- £3,  J3,  Pj,Rij  <  £' 


Using  the  scales  of  the  physical  parameters  derived  above,  we  present  the  optimized 
model  considering  only  leading  order  and  first  order  terms.  The  electromagnetic  model 
simplifies  to  2-D  wave  equation.  By  assuming  magnetic  walls  boundary  condition  and 
jcot  time  dependence,  the  electromagnetic  fields  are  approximated  using  a  cavity 
resonance  model.  The  optimized  model  is  given  by: 


d2ez  (PeT  2 

fo5~+ ez  =  JfajJ* 
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An  analytical  representation  of  the  solution  in  terms  of  eigen-function  expansion  is  given 
as: 


El  <  JH  lfm  n  > 
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mn 


As  inertia  effects  are  negligible  ,the  mechanical  model  is  reduced  to  linear  elastostatic 
problem.  The  effects  of  gravity,  electromagnetic  coupling  force  and  thrust  are  dominant. 
Deformation  along  z  is  dominant  due  to  thin  plate  approximation.  By  using  Hooke's  law 
as  the  constitutive  relation  and  applying  the  same  procedure  as  above,  a  leading  order 
approximation  for  determination  of  deformation  along  z  is  given  by: 
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where,  K  is  a  material  constant. 


It  should  be  noted  that  these  solutions  work  as  an  initial  estimate  to  the  solution.  For 
better  results  we  need  to  consider  second  order  corrections  as  well.  We  now  look  into  the 
behavior  of  a  regime  where  some  of  the  parameters  are  varied  by  order  of  s  and  study 
how  the  governing  equations  change  accordingly.  Consider  the  regime  where  order  of 
magnitude  of  electro-magnetic  fields  is  much  smaller  with  frequencies  of  the  order  of 
MHz  .  The  wave  equation  governing  electromagnetic  waves  is  now  modified  to: 
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Thermo-mechanical  problem:  We  have  used  second  law  of  thermodynamics  and 
Clausius-Duhem  inequality  to  perform  a  complete  combinatorial  analysis  of  72  possible 
energy  functions  to  describe  the  coupled  multi-physics  behavior.  Constitutive  equations 
have  been  derived  from  these  energy  potentials  in  a  rigorous,  systematic  fashion. 

Thermo-Electro-Magneto-Mechanical  (TEMM)  Constitutive  Equations: 

The  standard  independent  variable  pairs  and  the  corresponding  potential  functions  for  a 
thermo-mechanical  problem  are  listed  in  the  table  below: 


Independent  Variables 

Potential  function 

(F,n) 

£  -  Internal  energy 

(F,  0) 

i|i  -  Helmholtz  free  energy 

(p,  n) 

X-  Enthalpy 
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(P,  6) 


4>  -  Gibbs  free  energy 


The  Helmholtz  free  energy  is  a  function  of  mechanical  independent  variable  F  and 
thermal  independent  variable  .  There  are  9  such  potentials,  for  instance: 

tl>  =  e,  m) 

where,  the  independent  electrical  and  magnetic  variables  in  this  case  are  e  and  m, 
respectively.  Calculating  the  derivative: 
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This  equation  should  hold  for  all  processes,  which  gives: 
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The  remaining  dependent  variables  (electric  displacement  d  and  magnetic  flux  b)  are 
deduced  from  the  solution  via  post-processing.  Adherence  to  this  framework  ensures 
consistency  of  the  characterization  of  the  coupled  material. 

There  are  8  pairs  of  independent  variables  possible  for  a  thermo-mechanical  problem. 
With  three  electric  variables  (e,  d,  p)  and  three  magnetic  variables  (h,  b,  m),  9 
combinations  are  possible  for  each  of  these  8  pairs. 


In  all,  there  are  seventy-two  potentials  possible  for  TEMM  material 
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Task  2:  A  Semi-Analytical  Formulation  for  Multi-Scale  Element  Development  in  EM  Formulation: 

In  collaboration  with  Prof.  Ghosh  and  his  team  we  have  worked  on  development  of 
Multi-Scale-Super-Element  for  homogenization  of  microstructural  RYE’s  with 
parametric  forms  of  boundary  conditions  with  the  following  properties: 

Why  Super  Element? 


1)  Shape  functions  derived  from  analytical  solutions  increases  accuracy  of 
simulation 

2)  Inherently  accommodates  multi-scale  inhomogeneity 

3)  Fine  discretization  of  conducting  wire  is  robust,  does  not  take  physics  into 
consideration  whereas  super  element  method  provides  a  superior  modeling 
technique 

4)  Dependence  on  analytical  solutions  allows  for  parameter  study 


5)  Super  Element  exploits  shape  functions  derived  from  semi-analytical  solutions  in 
improving  accuracy  and  enabling  parameter  studies  for  design  and  optimization 
problems 

6)  Two-media  problem  is  reduced  to  single  medium  using  volume  equivalence  and 
surface  equivalence  principles,  thus  accommodating  multi-scale  inhomogeneity 


We  have  successfully  applied  the  super-element  to  a  simple  1-D  and  2-D  problems  and 
deduced  the  integral  equations  governing  multi- scale  EM  behavior  and  we  are  currently 
working  on  integrating  it  with  3-D  problem 
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Research  Summary 

The  UCLA  team  focused  on  research  to  enable  a  new  paradigm  of  structurally  integrable 
antenna  electronics.  The  concept  is  to  use  magneto-electric  materials  for  transmission  and 
reception  surfaces  (antennae)  of  broadcast  signals.  This  would  not  only  replace  existing 
conductive  materials  in  antenna  systems  but  could  allow  phasing  control  circuitry  to  be 
integrated  in  to  the  antenna  surfaces  of  an  antenna  array.  Magneto-electric  (ME) 
materials  exhibit  unique  behavior  in  a  number  of  aspects  due  to  an  inherent  cross 
coupling  of  electric  and  magnetic  fields  with  the  material  structure.  Due  to  this  cross 
coupling,  ME  materials  will  become  electrically  polarized  in  an  applied  magnetic  field 
and  conversely  will  exhibit  magnetization  when  exposed  to  an  electric  field. 

Among  the  interesting  phenomena  seen  in  ME  materials  are  strong  dependence  on 
propagation  direction  for  the  determination  of  a  waves  phase  velocity  and  its  modes  of 
propagation.  The  speed  modulation  capability  of  the  materials  may  prove  useful  in  the 
development  of  electrically  small  antennas.  Further,  since  the  phase  velocity  can  vary 
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significantly  for  each  propagation  mode,  polarization  rotation  commonly  occurs  in  ME 
materials,  providing  the  potential  for  polarization  modifying  devices.  Also,  due  to  the 
anisotropy  of  the  materials,  wave  polarization  can  be  measured  directly  with  an 
appropriate  electrode  configuration  as  opposed  to  the  need  to  filter  polarization 
components  and  then  measure  intensity.  Another  significant  factor  of  typical  ME 
materials  is  that  they  respond  to  various  field  biases.  For  example,  the  magnitude  of 
cross-coupling  parameters  changes  under  the  application  of  a  static  electric  field. 
Consequently,  ME  materials  are  considered  “active  materials”  which  can  be  tuned  on  the 
fly  by  application  of  an  appropriate  bias.  This  increases  the  possible  applications  of  the 
materials  and  may  prove  useful  for  phase  matching/altering  of  electrical  signals. 

The  applications  indicated  above  depend  on  the  ability  to  produce  ME  materials  with 
sufficient  cross-coupling.  Since  single  phase  ME  materials  generally  posses  weak  cross¬ 
coupling  parameters  their  use  at  this  time  is  impractical.  However,  strain  coupled 
composite  laminates  show  relatively  large  magneto-electric  coupling  effects.  In  addition 
to  providing  sufficient  coupling  characteristics,  laminate  construction  works  well  for  the 
thin  plate  and  thin  film  type  applications  envisioned  for  these  materials.  Therefore,  the 
work  conducted  under  this  program  focused  on  evaluating  the  coupling  effects  in  a  wide 
range  of  laminate  constructions  both  theoretically  and  experimentally. 

To  facilitate  use  of  ME  materials  in  the  previously  described  applications,  research 
efforts  at  UCLA  have  focused  on  three  main  areas:  Analytical  Wave  Propagation 
Modeling,  Development  of  the  Effective  (Homogenized)  ME  Material  Parameters  for 
Laminates,  and  the  Fabrication  Testing  and  Characterization  of  ME  structures  (including 
single  domain  entities).  Below  a  brief  synopsis  of  each  area  is  provided. 


Analytical  Wave  Propagation  Modeling: 

Wave  propagation  in  extended  media 

Our  work  has  focused  on  electro-magnetic  behavior  of  single-phase  homogenized 
structures  with  direct  cross-coupling  of  both  E  and  H  fields.  The  materials  analyzed  are 
assumed  to  be  lossless  although  the  formulations  used  for  calculations  work  for  general 
frequency  dependent  materials.  Utilizing  the  source  free  Maxwell’s  Equations 

Vxi  =  -B  V  x  H  =  -D  icy.  E  =  coB  k  x  H  =  -coD 

=> 

VB  =  0  V -3  =  0  k-B  =  0  k-D  =  0 

and  lossless  linear  magneto-electric  constitutive  relations  for  electric  flux  density  D  and 
magnetic  flux  density  B, 

E  =  kD  +  ~xB  and’  H  =  ~x  D  +  vB 

where  %  represents  the  second  order  magneto-electric  tensor.  Using  these  equations  we 
developed  a  generalized  wave  equation  that  allows  for  the  analytic  solution  of  plane  wave 
propagation  modes. 
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By  reducing  the  three  dimensional  problem  to  two  dimensions,  the  problem  is  made  more 
tractable  and  the  dispersion  relation  can  be  obtained  in  an  analytic  form.  Once  values  are 

CO 

determined  for  -  the  system  of  equations  can  be  solved  to  determine  the  modes  for  B. 

k 

The  B  solutions  then  determine  D  from  the  second  system  of  equations. 

The  dispersion  relation  is  a  4th  order  polynomial  in  the  inverse  of  the  index  of  refraction. 
For  special  cases,  i.e.  when  the  material  cross-coupling  goes  to  zero,  the  dispersion 
relation  can  be  reduced  to  an  effective  2nd  order  equation  which  is  representative  of 
conventional  material  systems.  However,  for  general  ME  material,  four  distinct  roots,  i.e. 
wave  phase  velocities,  can  exist  for  a  given  propagation  vector.  Generally  this  translates 
to  two  modes  of  propagation  in  the  forward  direction  and  two  different  modes  of 
propagation  in  the  opposite  direction,  which  is  quite  different  from  the  behavior  of 
conventional  materials.  That  is  the  material  is  non-reciprocal  and  forward  velocity 
propagation  speeds  are  distinctly  different  from  backward  velocity  propagation  speeds. 

For  the  case  of  a  “homogeneous”  magneto  electric  material,  we  have  examined  the  phase 
velocity  profiles  for  the  archetypal  cases.  Figure  1  shows  plots  of  phase  velocity  with 
angle  of  propagation  through  the  material  indicating  asymmetry  between  forward  and 
backward  travel,  a  characteristic  of  the  non-reciprocal  ME  materials.  For  each 
propagation  direction  there  are  two  modes  of  propagation  found,  as  exhibited  by  the  two 
phase  velocity  surfaces  shown.  For  a  general  propagation  vector  it  can  be  seen  that  the 
phase  velocity  of  each  mode  is  different.  This  velocity  difference  causes  a  polarization 
rotation  of  the  wave  as  it  propagates  through  an  ME  material.  The  forward  and  reverse 
rotations  are  different  in  general  as  can  be  seen  from  the  asymmetry  in  each  plot.  At  some 
corresponding  points  on  the  surfaces  (i.e.  at  the  same  angles),  we  see  the  velocities  are 
the  same.  At  these  directions  the  waves  propagate  without  rotation  and  the  incoming 
mode  is  unaltered. 
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Figure  1  Phase  velocity  profiles  as  a  function  of  wavevector  direction  for  a  material  with  Xu 
coupling.  Each  surface  represents  the  phase  velocity  of  one  propagation  mode. 

The  previous  plot  showed  a  basic  illustration  of  the  non-reciprocal  nature  of 
magnetoelectric  laminates.  However,  more  interesting  situations  have  been  observed  in 
the  parameter  space,  such  as  zero  velocity  modes  like  the  one  shown  in  figure  2.  For  large 
cross-coupling  parameters,  the  phase  velocity  in  a  certain  direction  will  approach  zero. 
Thermodynamic  limits  prevent  the  cross-coupling  from  actually  achieving  the  zero  point; 
however,  it  should  be  theoretically  possible  to  get  infinitesimally  close.  For  real  materials 
it  is  uncertain  how  close  but  large  velocity  ratios  should  be  realistic. 


Figure  2  Phase  velocity  profile  for  a  material  with  %l2  and  coupling.  The  point  structure 
shown  in  the  plot  indicates  a  region  in  which  the  phase  velocity  decreases  to  zero.  The  zero 
velocity  point  corresponds  to  a  theoretical  limit  but  arbitrarily  small  velocities  are  theoretically 
possible. 
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Our  work  shows  that  phase  shifting  and  polarization  rotation  applications  are  possible. 
The  separation  of  polarization  components  is  more  direct  with  magneto-electric  materials 
than  with  “conventional”  materials  since  field  polarization  components  may  be  separated 
by  using  material  coupling  rather  than  filtering  a  specific  polarization  and  then  measuring 
the  component  amplitude. 

Interface  behavior  of  ME  materials 

The  previous  work  was  extended  to  examine  the  behavior  of  waves  at  the  interface  of  two 
ME  materials  or  the  interface  between  an  ME  and  a  conventional  material.  While  Snell’s 
law  still  holds  for  ME  materials,  the  fact  that  the  wave  vector  is  dependent  on 
transmission  angles  through  the  material  creates  additional  interesting  behavior  not 
normally  observed.  The  propagating  EM  wave  is  split  on  reflection  (ME-ME  interface 
only)  and  transmission,  although  for  some  materials  parameters  there  are  propagation 
directions  which  do  not  result  in  bireffingent  behavior  on  transmission  (See  Figure  3). 
Our  current  work  has  been  to  develop  a  model  that  calculates  the  wave  vectors, 
reflection/transmission  coefficients,  and  field  components  of  an  incident  plane  wave  onto 
an  ME  interface  with  an  emphasis  on  identifying  interesting  “propagation”  modes  that 
either  rapidly  disperse  wave  energy  or  allow  forward  propagation  without  significant 
modification. 


Air  X  ME1  ME2 


Figure  3  Typical  Transmission/Reflection  characteristic  at  ME  interface. 

Wave  behavior  in  constituent  phases  of  ME  laminates 

To  better  understand  the  direct  elastic  interaction  effects  of  the  proposed  ME  laminates 
the  wave  equation  and  techniques  used  above  were  applied  to  the  piezoelectric  and 
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magnetoelectric  phases  individually.  Using  the  mechanical  equilibrium  equation  and 
small  strain  approximation  in  spatial  and  time  frequency  domain: 

k2  T 

e= - vNTNq, 

pco  =  = 

where  epsilon  is  the  strain,  sigma  is  the  stress  (both  in  Voigt  notation),  k  is  the  magnitude 
of  the  wavevector,  omega  the  angular  frequency  and  N  relates  to  the  direction  cosines  of 
the  unit  vector  in  the  direction  of  wave  propagation, 

nx  0  0  rc3  0  n2 

N  =  0  n2  0  0  o3  nx  with  the  ni  the  direction  cosines  of  the  wavevector  to 

0  0  n3  nx  n2  0 

coordinate  axis.  Using  the  piezoelectric  constitutive  equations  as  an  example, 
s  =  sq  +  dT  E,  D  =  da+  eE ,  and  B  =  pH_  ,  we  combine  the  equilibrium  equation  above 

with  the  first  constitutive  equation  to  obtain: 

(  k2  T  }  T 

I  - TNT  N-s\q  =  dT  E. 

\  pco  =  =  =) 

k 

The  matrix  on  the  left  hand  side  is  invertible  exactly  when  the  ratio  of  —  is  not  the  speed 

co 

of  a  mechanical  wave,  which  we  would  have  if  the  d  matrix  was  zero.  This  can  now  be 
used  to  eliminate  <r  from  the  other  constitutive  equation  to  form: 

(  k2  V1" 

D  =  eE  with  e  =  e  +  d  I  - t-NtN-s  I  dT  . 

=  =  =  =\pco  =  =  =)  = 

This  effectively  converts  the  electromagnetic -mechanical  problem  to  a  strictly 
electromagnetic  problem  which  can  be  handled  as  a  normal  eigenvalue  problem.  The 
piezomagnetic  case  can  be  handled  directly  by  duality,  switching  the  roles  of  p  and  e , 

etc. 

Once  the  field  quantities  are  solved  in  their  respective  piezoelectric  and  piezomagnetic 
materials,  boundary  conditions  can  be  applied  at  an  interface  to  solve  for  any  propagating 
wave  modes.  Solutions  for  propagating  waves  are  analytic  or  at  least  soluble  to  arbitrary 
precision  and  can  also  be  used  to  test  the  homogenization  models  for  ME  laminates 
proposed  below. 
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Development  of  effective  laminate  properties: 


A  survey  of  available  ferroelectric  and  ferromagnetic  materials  was  performed  to  identify 
optimal  combinations  to  maximize  the  magnetoelectric  coupling.  The  primary  purpose  of 
the  survey  was  to  identify  ferroelectric  material  with  large  electromechanical  coupling 
and  ferromagnetic  material  with  large  magnetomechanical  coupling  as  well  as  the 
permittivity  and  permeability  of  the  respective  materials.  Following  this  survey,  an 
analytic  method  was  developed  to  determine  the  effective  material  properties  of 
magnetoelectric  laminate  structures. 

The  constitutive  equations  of  any  given  layer  can  be  expressed  in  the  following  form. 
Where  y  is  strain,  D  is  electric  flux,  B  is  magnetic  flux,  a  is  stress,  E  is  electric  field,  H  is 
magnetic  file,  S  is  the  compliance  matrix,  d  is  the  piezoelectric  tensor,  q  is  the 
piezomagnetic  tensor,  s  is  permittivity,  p  is  permeability,  and  m  is  the  magnetoelectric 
tensor. 


Rearranging  the  equations  allows  us  to  relate  stress,  electric  flux,  and  magnetic  flux  to 
strain,  electric  field  and  magnetic  field  in  each  layer. 
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By  integrating  through  the  thickness  of  the  layered  structure  (z)  an  expression  is  found  to 
relate  mid-plane  strain  (y°),  curvature  (k),  electric  filed  and  magnetic  field  of  the  laminate 
to  the  resultant  forces  (N),  resultant  moments  (M),  electric  flux  and  magnetic  flux. 
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The  matrix  in  the  above  relation  represents  the  effective  material  properties  for  the 
laminate  as  a  whole.  These  effective  laminate  properties  are  used  in  Maxwell’s 
Equations  (previous  sections)  to  determine  EM  wave  propagation  in  the  laminate. 
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There  is  a  problem  with  the  above  formulation  in  that  it  essentially  ignores  out  of  plane 
strain  components  which  may  be  important,  particularly  for  wave  propagation 
applications.  UCLA  is  working  on  an  alternative  formulation  which  may  better 
approximate  the  homogenized  material  parameters  by  taking  into  account  the  out  of  plane 
interactions  that  arise  from  piezoelectric  and  piezomagnetic  coupling. 

If  we  take  the  z  axis  as  the  out  of  plane  for  a  laminate  ME  material  we  find  that  half  of 
the  field  components  are  continuous  across  all  layers  of  the  composite  as  required  by 
boundary  conditions  of  the  respective  mechanical  or  electromagnetic  governing 
equations.  For  the  current  situation  this  implies 

( yl,y2,cr3,a4,a5,y6,El,E2,D3,Hl,H2,B3 )  are  continuous  quantities  across  all  interfaces 
while  (crl,cr2,y3,y4,y5,cr6,Dl,D2,E3,Bl,B2,H3)  generally  change  discontinuously  across 
the  interfaces. 

For  linear  materials  it  is  straightforward  to  rearrange  the 
{a  D  B  J=[S]{  y  E  B  constitutive  form  shown  above  to  a  hybrid  form 
such  that : 


{crl,(T2,y3,Y4,r5,o-6,Dl,D2,E3,Bl,B2,H3}T  =  \M]{Yl,r2,cr3,(J4,(J5,y6,El,E2,D3,HvH2,B3}T 


where  [M]  represents  the  rearranged  constitutive  parameter  matrix.  For  a  layered 
material  we  have  constitutive  equations  of  the  above  form  for  each  layer  where  [M] 
contains  the  material  specific  dependencies  for  each  layer.  If  we  denote  the  above 
equation  for  the  nth  layer  as 

Y(n)  =  [M](n)^n) 

where 

Y(n)  =  7l,cr2,y3,y4,y5,a6,Dl,D2,E3,Bl,B2,H3  }r  for  the  nth  layer  and 

^n)  =  {ri,r2,°3,V4’0'5’r6’El’E2’D3’Hl’H2’B3}T 

then  the  integral  through  the  thickness  of  the  material  becomes 

h+  h * 

|  dzY(n)  =  jdz[Mp^n). 

h~  h~ 

Assuming  each  layer  is  homogeneous,  dividing  by  the  thickness  of  the  composite 
t  =  h+  —  h~  gives 
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<r)=;Z[Mf>  j  dzX^n) . 

^  n= 0  z« 

The  angle  brackets  indicate  the  average  in  z  of  the  field  quantities  contained  in  Y .  We 


note  that  for  a  1st  order  estimate  where  each  field  quantity  is  constant  for  each  layer,  the 
above  simplifies  to 


<r>= 


N  \ 

X  ,  vf  =  volume  fraction  of  the  nth  layer. 

vn=0  ' 


Since  the  field  variables  on  the  right  side  are  continuous  they  must  be  constant 
throughout  the  laminate  thickness  to  first  order.  If  we  then  rearrange  the  resulting  matrix 

equation  into  |  <j  D  B  J  =  [S']  y  E  B  format,  [5]  being  a  fairly  involved 

mixture  of  constitutive  parameters  from  each  layer,  we  reduce  to  the  cross-coupling 
estimations  seen  commonly  in  the  literature  for  bi-  or  tri- layer  structures  when  N=2  or  3. 
The  form  suggested  above  provides  a  convenient  manner  for  calculating  the  cross¬ 
coupling  for  multi-layer  systems  compared  with  the  approach  conventionally  employed. 

Additionally,  for  higher  order  estimates,  the  field  quantities  contained  on  the  right  side  of 
the  integral  equation  can  be  estimated  by  series  expansions  compatible  with  the 
continuity  of  the  underlying  fields.  In  particular,  we  should  be  able  to  express 
X  =  X(Kx  (x,y),  K2  (x,y),z)  where  Kx  and  K2  represent  the  local  curvatures  at  a 
reference  plane  in  the  laminate  and  express  the  constitutive  relations  for  the  plate  as  a 
function  of  the  resultant  stresses,  moments,  curvature,  and  averaged  field  quantities  of  the 
composite. 

Fabrication,  Characterization  and  Testing  of  ME  materials 


Comparison  of  Effective  Direct  and  Converse  Magneto-electric  Effects  in  Laminate 
Composites 

There  is  a  good  deal  of  confusion  in  the  general  literature  describing  direct  and  converse 
magneto-electric  effects.  This  stems  largely  from  a  lack  of  standardization  in  what 
physical  quantities  are  reported  and  the  specific  methods  of  measurement  (ie  -  whether 
the  entire  laminate  response  is  reported  or  only  the  response  in  say  the  piezomagnetic 
phase).  Theoretically  from  thermodynamic  considerations  one  should  expect  that  the 

dD 

direct  magneto-electric  effect  (DME)  defined  as  —  and  the  converse  magneto-electric 

dH 

dB 

effect  (CME)  defined  by  — should  be  equivalent.  An  experiment  was  conducted  on  a 

dE 

tri-layer  PZT-TerfenolD-PZT  plate  to  compare  actual  DME  and  CME  coefficients  with 
expected  results  from  an  analytic  model  based  on  the  constituent  phases  indicated  above. 

Our  analysis  on  the  tri-layer  structure  predicts  a  CME  coefficient,  ac ,  and  a  DME 
coefficient,  ac ,  as  follows; 
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aD  = 


dD3 

dH \ 

Sl  Jm  +  r^Sl\tp 

■tm  and  ac  = 


dBl 

^31, pd\\,m 

dE3 

e E  f  t 

**1  Vm  ^  A^\\Lp 

•2 1 


D3  and  E3  are  the  electric  displacement  and  electric  field  along  the  piezoelectric  thickness 
direction,  respectively.  Bi  and  Hi  are  magnetic  flux  density  and  magnetic  field  strength 
along  the  length  direction,  respectively,  .vf,  and  ,s^  are  the  compliance  coefficients  at 
constant  electric  field  and  constant  magnetic  field,  respectively.  d3lp  and  dn  are  the 
piezoelectric  and  piezomagnetic  coefficients,  respectively. 

To  evaluate  the  theoretical  analysis,  a  trilayer  magnetoelectric  laminate  was  fabricated. 
The  laminate  consisted  of  a  single  layer  of  (112)-oriented  Tbo.3Dyo.7Fe1. 92  (Terfenol-D) 
magnetostrictive  alloy  plate  (i.e.,  oriented  along  the  1 -direction)  between  two  layers  of 
oppositely  poled  Pb(ZrxTii.x)03(PZT-5H)  piezoelectric  plates.  The  (112)-oriented 
Terfenol-D  plate  was  obtained  from  ETREMA  Inc.  (Ames,  IA),  and  had  a  dimension  of 
24.45x6.7x1. 6mm  .  The  piezomagnetic  coefficient,  dU  m ,  of  the  (112) -oriented  Terfenol- 

D  plate  was  measured  to  be  5-7  nm/A  under  optimum  magnetic  bias  field  and  the  elastic 
compliance  coefficient  was  30-40  pm2/N.  The  PZT-5H  plates  were  obtained  from  Piezo 
Systems,  Inc.  (Cambridge,  MA),  and  had  a  dimension  of  24. 5x6. 7x0. 267mm3.  The 
piezoelectric  coefficient,  d31  ,  and  elastic  compliance  coefficient,  .vf, ,  of  the  PZT-5H 

plate  are  -320pC/N  and  16.5  pm2/N,  respectively.  The  Terfenol-D  plate  was  sandwiched 
between  PZT-5H  plates  with  a  7 6 /.an  thick  silver  conductive  film  (CF3350,  Emerson  & 
Cuming).  The  laminate  was  cured  at  125  C  for  2  h  with  mechanical  load  to  enhance  the 
bonding. 

The  DME  effect  of  the  laminate  sample  was  characterized  by  measuring  magnetic-field- 
induced  electric  polarization  (MIE).  An  electromagnet  was  employed  to  provide  a  DC 
magnetic  field  bias  ranging  from  0  to  3000  Oe.  A  small  custom  made  solenoid  coil 
applies  a  small  AC  magnetic  field  (0.1  to  2  Oe).  The  induced  out  of  plane  ME  voltage, 
through  the  thickness  of  the  PZT  layers,  was  measured.  The  CME  effect  was 
characterized  by  measuring  electric-field-induced  magnetization  (EIM).  A  sine  wave 
reference  source  coupled  with  a  voltage  amplifier  applied  an  AC  electric  voltage,  from  1 
to  100  V,  across  the  PZT  layers.  The  induced  magnetization  in  the  Terfenol-D  layer  is 
calculated  from  the  emf  voltage  induced  in  a  search  coil,  which  was  measured  with  a 
lock-in  amplifier.  It  is  important  to  note  that  the  search  coil  was  wrapped  around  the 
entire  laminate  and  not  just  the  magnetostrictive  layer.  This  approach  is  typically  used  in 
reporting  CME  response. 

We  found  the  effective  DME  and  CME  coefficients  to  be  of  the  same  order  of  magnitude, 
but  were  significantly  different  (see  Figure  4).  The  measured  values  were  also 
significantly  smaller  than  predicted  by  analysis. 


71 


12 

10 

8 

6 

4 

2 

0 


Figure  4  Effective  DME  and  CME  coefficients  and  the  ratio  of  two  effective  coefficients 
as  a  function  of  excitation  frequency.  The  inset  shows  details  of  low-frequency  data. 


The  experimental  results  suggest  that  closer  scrutiny  of  these  results  and  other  data  in  the 
literature  is  needed  to  better  report  measurement  values  that  can  be  compared. 


Influence  of  electric  voltage  bias  on  converse  magnetoelectric  coefficient  in 
piezofiber/Metglas  bilayer  laminate  composites 

Magneto-electrics  are  active  materials  whose  behavior  is  strongly  influenced  by  bias 
fields  applied  to  the  materials  under  test.  For  example  magnetic  fields,  electric  fields  and 
stress/strain  fields  all  can  be  used  to  control/adjust  ME  material  properties.  However 
results  in  the  literature  have  only  focused  on  the  influence  of  magnetic  fields  and  there  is 
a  virtual  absence  of  experimental  data  on  the  influence  of  either  electric  fields  or 
stress/strain  bias  fields.  Therefore,  experiments  were  conducted  on  a  piezofiber/Metglas 
bilayer  (PFM)  composite  subjected  to  both  a  dc  electric  voltage  bias  and  subsequently  a 
mechanical  bias. 


The  tests  done  here  focused  on  evaluating  the  optimum  electric  bias  fields  on  a  PFM 
fabricated  at  UCLA.  The  piezofiber  layer  is  a  macrofiber  composite  (MFC)  from  Smart 
Material  Corp.,  with  dimensions  of  28x7x0.3  mm3.  The  piezoelectric  constant  J33  value  is 
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Ratio  of  two  effective  coefficients 


approximately  0.75  ppm/V  and  the  Young’s  modulus  is  30  GPa  in  the  fiber  direction.  The 
Metglas  layer  is  a  2826MB  type  iron  nickel  based  alloy  (Metglas  Inc.),  with  dimensions 

o 

of  28x7x0.025  mm  .  The  PFM  layers  are  bonded  together  with  Loctite®  Epoxy.  The 
measured  thickness  of  the  entire  sample  is  approximately  0.375  mm. 

To  characterize  the  CME  response,  a  dc  magnetic  field  bias  7/bias,  up  to  40  Oe,  is  applied 
along  the  longitudinal  direction  of  the  sample  with  a  Helmholtz  coil.  A  separate  dc 
electric  voltage  bias  Vbias  ,  -300  to  700  V,  is  applied  with  a  dc  source  coupled  with  a 
voltage  amplifier.  An  ac  electric  voltage  Vac  ,  0.1-60  VrmS,  is  applied  using  a  sine  wave 
generator  coupled  to  the  voltage  amplifier.  The  induced  magnetic  field  is  measured  with  a 
coil  wrapped  around  the  sample. 

Previous  studies  on  similar  materials  show  that  a  dc  magnetic  field  bias  can  be  used  to 
maximize  the  CME  coefficient.  Our  experimental  data  also  show  that  the  optimum  dc 
magnetic  field  bias  is  a  function  of  applied  dc  electric  voltage  (see  Figure  5). 
Furthermore,  an  optimum  dc  electric  voltage  bias  exists  which  maximizes  the  converse 
magnetoelectric  coefficient. 
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Figure  5  Measured  magnetic  flux  density  B  as  a  function  of  magnetic  field  bias  FI  bias  with 
different  electric  voltage  biases. 


Our  data  indicates  that  the  CME  effect  is  strongly  dependent  on  dc  electric  voltage  bias 
with  the  value  changing  by  up  to  70%  over  the  dc  electric  voltage  bias  values  studied. 


Influence  of  mechanical  load  bias  on  CME  effect  in  ME  laminates 


UCLA’s  previous  studies  have  focused  primarily  on  the  bias  effects  of  electric  and 
magnetic  fields  on  the  performance  of  ME  materials.  Of  equal  or  greater  importance  is 
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understanding  the  influence  of  mechanical  loading  on  ME  laminates.  In  real  world 
settings,  ME  laminates  will  undergo  mechanical  biasing  both  from  intended  application 
of  loads  as  well  as  in  response  to  environmental  factors.  Consequently,  data  is  needed  to 
characterize  ME  materials  under  diverse  loading  conditions. 

Here  a  piezofiber/Metglas  (PFM)  ME  laminate  has  been  integrated  into  a  graphite  epoxy 
composite  (GEC)  to  study  the  CME  effect.  Experimental  data  on  a  PFM/GEC  driven  with 
an  ac  electric  voltage  subjected  to  both  a  dc  magnetic  field  bias  and  a  dc  mechanical  load 
bias  are  presented.  Results  of  these  tests  indicate  that  both  the  mechanical  load  and  the  dc 
magnetic  field  strongly  influence  the  CME  response.  Furthermore,  an  optimum 
mechanical  load  exists  to  maximize  the  CME  coefficient.  These  results  reveal  that  the 
CME  coefficient  can  be  further  increased  with  a  proper  mechanical  load  bias.  Therefore, 
the  selection  of  an  appropriate  mechanical  preload  as  well  as  dc  magnetic  bias  will 
maximize  the  CME  response  and  sensitivity  in  ME  laminates  as  well  as  integrated 
structural  systems. 

A  PFM/GEC  laminate  composite  was  fabricated  for  the  tests  performed.  The  piezofiber 
layer  is  a  macrofiber  composite  (Smart  Material  Corp.)  with  dimensions  of  28x7x0.3 
mm  .  The  piezoelectric  constant  r/33  value  is  approximately  0.75  ppm/V,  and  Young’s 
modulus  is  28  GPa  in  the  fiber  direction.  The  Metglas  layer  is  a  2826MB  type  iron  nickel 
based  alloy  (Metglas  Inc.),  with  dimensions  of  28x7x0.025  mm  .  The  Young’s  modulus 
of  Metglas  is  110  GPa.  The  graphite  epoxy  AS4/  3501-6  composite  (GEC)  is  a  [(0°  /+- 
45°  /90°)5]3  quasi-isotropic  layup  with  an  effective  in-plane  Young’s  modulus  of  50  GPa. 
The  GEC  dimensions  are  160x49x7  mm3.  The  PFM  ME  laminate  is  attached  to  the  center 
of  one  side  of  the  GEC  using  Loctite®  epoxy.  For  analytical  modeling,  the  in-plane 
Poisson’s  ratio  for  the  MFC  is  0.36  and  for  GEC  layer  is  0.3. 

The  PFM/GEC  composite  sample  was  mounted  in  a  load  frame,  and  both  tensile  and 
compressive  static  loads  applied.  The  PFM  bilayer  is  approximately  40  cm  away  from  the 
gripped  region  of  the  GEC.  Two  strain  gauges  are  attached  on  both  sides  of  the 
PFM/GEC  sample.  To  characterize  the  CME  response,  an  ac  electric  driving  voltage  Vac 
is  applied  on  the  MFC  layer  using  a  sine  wave  source  fed  to  a  voltage  amplifier.  A  dc 
magnetic  field  bias,  up  to  40  Oe,  is  applied  along  the  longitudinal  direction  of  the  sample 
using  a  Helmholtz  coil.  The  induced  magnetization  in  the  PFM/GEC  sample  is  measured 
using  a  120-turn  search  coil  wrapped  around  the  entire  PFM/GEC  composite. 

A  plot  of  obtained  magnetization  data  at  seven  different  mechanical  loads/strain  values  is 
shown  in  Figure  6.  Our  experimental  results  indicate  that  the  mechanical  load  bias 
strongly  influences  the  value  of  the  optimum  //bias-  A  proper  preload  should  be  selected  to 
optimize  the  CME  response  and  sensitivity  for  magnetic  field  detection  in  both  ME 
laminates  as  well  as  integrated  structural  systems. 
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Figure  6  Normalized  electric -field-induced  magnetization  ( CME  effect )  as  a  function  of 
dc  magnetic  field  bias  under  different  mechanical  loads/strains. 


Nanostructure  Fabrication 

To  better  understand  the  response  of  magnetoelectric  in  the  small  scale  and  thin  film 
structure,  we  fabricated  an  ME  thin  film  material.  As  previously  reported  the  deposition 
of  both  Ni  on  PZT  thin  film  were  on  silicon  substrates  with  the  Ni  film  is  approximately 
100  nm  thick  and  the  PZT  film  is  about  1.28  micron.  A  number  of  tests  were  conducted 
on  the  structure  below  we  describe  the  magnetic  force  microscope  MFM  results  when  a 
voltage  is  applied  to  the  PZT.  During  this  test  the  surface  of  the  ME  thin  film  was 
monitored  with  an  MFM  and  a  typical  stripe  domain  pattern  on  the  Ni  phase  of  the 
bilayer  film  appeared  without  an  applied  voltage.  After  application  of  10V  to  the  PZT 
phase  of  the  structure  MFM  imaging  reveals  significant  differences  in  the  stripe  domain 
patterns  (see  Figure  7  below).  These  images  confirm  significant  ME  effects  in  nano- 
structured  materials  using  strain  coupled  piezoelectric  and  piezomagnetic  phases.  When 
an  electric  field  is  applied  to  the  PZT-film,  mechanical  coupling  between  the  PZT-film 
and  the  Ni-film  induces  strain  in  the  Ni.  The  magnetostrictive  properties  of  the  Ni  then 
induce  changes  in  the  observed  domain  pattern.  The 
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Figure  7.  Magneto-electric  effect  in  a  Ni/PZT  thin  film  material.  Image  at  left  shows 
sample  with  no  voltage  applied.  The  image  at  right  is  the  same  sample  with  10V  ( approx . 
1  MV/m  field).  Changes  in  stripe  domain  pattern  is  evident  by  comparison  of  labeled 
regions. 


Next  we  fabricated  a  magnetoelectric  nanostructure  and  demonstrated  that  electric  control 
of  the  magnetic  spin  structure  is  possible.  An  illustration  of  the  magnetoelectric 
nanostructure,  a  layered  Ni  nanobar  on  a  PZT  thin  film  structure,  is  shown  in  Figure  8. 

The  MFM  images  indicate  substantial  change  in  the  magnetization  of  the  nanobar  with 
the  application  of  voltage  to  the  piezoelectric  phase  of  the  ME  structure.  The  length  and 
width  of  the  Ni  nanobar  is  approximately  380  nm  xl50  nm.  As  in  the  previous  structure, 
strain  coupling  between  phases  produces  the  changes  in  magnetization  seen  in  the  Ni. 


Figure  8.  Illustration  of  change  in  magnetization  structure  seen  with  MFM.  The  image 
on  the  left  is  with  zero  volts  and  the  image  on  the  right  shows  1.5  V  applied  to  the  PZT 
phase  of  the  ME  nanostructure. 
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Our  tasks  for  the  GameChanger  project  were  in  two  areas,  one  short  term  that  concerns 
the  use  of  magnetic  nanowires  in  antenna  systems  and  one  long  term  that  seeks  to 
develop  general  methods  to  calculate  the  properties  of  magnetic  and  ferroelectric 
properties  of  nanostructures  from 

first  principles.  The  former  concerns  the  use  of  nanowires  to  increase  Sr  and  pr 
proportionately  to  slow  the  wave  while  matching  impedance.  The  key  point  about  the  use 
of  nanowires  is  that  they  can  exist  in  the  single  domain  state,  without  domain  walls.  The 


presence  of  domain  walls  has  previously  been  linked  to  loss  in  ferromagnetic  elements. 
Thus,  it  is  critical  to  make  predictions  of  the  critical  radius  of  a  nanowire  at  which  the 
single  domain  state  loses  stability,  and  especially  to  understand  how  this  depends  on  the 
material  constants.  This  will  guide  the  selection  of  the  best  kind  of  magnetic  material  for 
this  possible  application.  The  second  thrust  of  our  research  is  a 

long  term  effort  to  develop  a  first  principles  method  for  the  systematic  calculation  of 
properties  of  nanostructures,  and  to  use  this  to  seek  new  nanostructures  with  desirable 
properties  for  antennas. 

On  the  first  area  of  research,  we  have  developed  general  methods  for  calculating  the 
breakdown  of  the  single  domain  state  in  nanowires.  It  is  based  on  a  mathematical  analysis 
of  the  full  energy  of  micromagnetics  including  exchange,  anisotropy,  demagetization, 
magnetostrictive  and  applied  field  contributions  to  the  energy.  For  this  we  have  given 
upper  and  lower  bounds,  and  we  have  also  developed  a  code  for  direct  calculation  of  the 
eigenfunction  at  breakdown. 


Typical  calculated  values  are  shown  in  Figure  1.  Typical  diameters  for  breakdown  of  the 
single  domain  state  in  nanowires  are  25-150  nm  with  little  dependence  on  aspect  ratio  for 
wire-like  structures. 
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Fig.  1  Rigorous  upper  and  lower  bounds  for  radius  at  breakdown  of  the  single  domain 
state. 

We  noticed  that  many  nanostructures  are  “objective  structures”,  that  is,  they  are 
composed  of  molecules  such  that  corresponding  atoms  in  each  molecule  see  the  same 
environment  up  to  orthogonal  transformation.  For  example,  a  single- walled  carbon 
nanotube  (of  any  chirality)  is  such  a  structure  with  one  atom  per  molecule. 
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During  the  past  year  we  have  made  major  progress  on  these  long  term  issues.  Our  main 
advances  have  been: 

1.  Major  progress  on  a  DFT  strategy  for  objective  structures,  including  advances  of  the 
basic  theory.  Our  methods  have  concentrated  on  the  development  of  an  analog  of  the 
Kohn-Sham  iterations.  For  objective  structures  such  as  buckyballs  and  rings  of  atoms, 
we  have  developed  a  real- space  finite  difference  algorithm  that  computes  the  electronic 
structure  while  taking  into  account  the  symmetry  that  underlies  these  structures.  For 
objective  structures  such  as  helical  rods  and  nanotubes,  we  have  decided  not  to  design  a 
real  space  strategy,  but  rather  to  develop  a  method  that  more  precisely  uses  the  isometry 
groups  associated  with  these  objective  structures  through  a  spectral  method.  We  have 
developed  an  analog  of  the  Bloch  theorem  for  objective  structures  for  some  of  the  groups 
describing  objective  structures.  We  have  also  developed  an  efficient  method  for 
generating  fundamental  domains  and  discretizing  them  as  well  as  a  pseudo -spectral  code 
that  has  been  used  for  producing  some  results.  A  description  of  this  work  will  be  given  at 
the  review. 

2.  We  have  developed  a  method  of  objective  molecular  dynamics  (OMD).  This  is  an 
exact  method  of  molecular  dynamics  —  every  atom  (typically  infinitely  many)  satisfies 
exactly  F=ma  —  but  only  a  few  atoms  need  be  simulated.  (The  time-dependent  positions 
of  other  atoms  are  given  by  an  explicit  formula.)  Certain  kinds  of  macroscopic  motions 
can  be  simulated.  As  an  application  of  the  method  we  have  simulated  the  behavior  of  a 
(6,6)  carbon  nanotube  pulled  at  constant  macroscopic  strain  rate,  one  of  the  possibilities 
allowed  by  OMD.  For  the  atomic  forces  we  have  used  Tersoff  potentials  for  carbon.  We 
have  been  able  to  reach  strain  rates  that  are  experimentally  accessible,  a  rarity  in 
molecular  dynamics  simulations.  From  a  series  of  50  simulations  at  various  strain  rates, 
we  have  obtained  the  following  unexpected  results: 

Several  modes  of  failure  were  observed,  including  fibrous  fracture  (the  nanotube  splits 
along  its  length),  cavitation  (holes  open  up  in  the  side  of  the  nanotube),  cross-sectional 
collapse  (the  tube  collapses  into  two  parallel  graphene  strips),  but  in  no  cases  have  we 
observed  the  expected  mode  of  failure  of  the  propagation  of  Stone-Wales  defects. 
Increasing  strain  rate  and  initial  temperature  paradoxically  led  to  lower  ductility,  Figure 
2.  In  those  cases  exhibiting  failure  at  low  strains  (~1%),  failure  was  preceded  by  large 
amplitude  vibrations  of  the  cross-section  at  a  frequency  of  about  100GHz. 
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Fig  2.  Strain  at  failure  of  (6,6)  carbon  nanotubes  as  a  function  of  strain  rate  and 
temperature  from  38  simulations 


Examples  of  these  large  amplitude  vibrations  of  the  cross-section  are  seen  in  Figure  3  in 
two  snapshots  from  a  simulation  at  high  strain  rate  and  initial  temperature  of  500K. 
This  also  shows  what  could  be  called  Stone-Wales  partials,  which  however  quickly 
healed  and  did  not  seem  to  contribute  to  plasticity  or  failure. 


Fig.  3.  Large  amplitude  vibrations  of  the  cross-section  of  a  carbon  nanotube  pulled  at 
constant  strain  rate  and  intitial  temperature  of  500K.  The  tan  atoms  show  a  loss  of 
coordination  and  the  formation  of  a  7-5  defect  (Stone- Wales  partial)  which  healed  in 
lOps. 
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Some  results  from  our  Objective  Density  Functional  Theory  codes  are  presented  below. 
Fig  4  shows  the  electronic  density  of  a  structure  which  has  reflection  symmetries  about 
the  3  principal  planes.  A  6th  order  differencing  scheme,  with  a  grid  size  of  0.2  a.u.  and 
603  gridpoints  was  used  to  solve  the  problem.  Results  (ground  state  energies,  etc) 
compare  well  with  KSSOLV  (Planewave  DFT  Code  developed  at  LBNL).  For  similar 
convergence  parameters,  CPU  Times  used  by  our  code  are  significantly  lesser  compared 
to  standard  finite  difference  and  plane-wave  DFT  codes. 


Fig  4:  “Proof  Of  Concept”  computation  of  electronic  structure  using  the  Objective 
DFT  algorithm:  Ring  of  8  aluminum  atoms  using  Heine  Abarenkov 
Psedupotential  for  Aluminum.  The  objective  domain  and  two  iso-surface  levels  of 
electronic  density  have  been  plotted. 


Some  results  for  the  Objective  Density  Functional  Theory  algorithm  for  helical  rods  and 
nanotubes  are  also  shown  below  in  Fig  5.  The  results  are  from  a  pseudo-spectral  code  that 
was  used  for  integrating  the  governing  equations.  We  have  develoved  techniques  for 
automatically  generating  fundamental  domains  for  these  objective  structures  and  our 
recent  success  has  been  the  development  of  a  spectral  method  that  systematically  exploits 
the  underlying  symmetries  of  these  structures. 


Fig  5:  Left  and  Middle  Images:  Electronic  charge  density  (two  levels  of  iso-surfaces)  and 
contours  plotted  over  two  fundamental  domains  for  an  Aluminum  nanotube  with  a  5 
degree  angle  of  twist.  Right  Image:  Charge  density  plot  across  the  cross  section  of  the 
same  tube  shows  a  large  delocalized  body  of  electrons.  The  results  are  from  a  pseudo - 
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spectral  code  that  was  used  to  solve  Kohn-Sham  equations  with  Bloch-like  boundary 
conditions. 

Our  theoretical  work  on  magnetic  materials  for  the  GameChanger  project  has  suggested 
an  interesting  new  material  that  achieves  unusually  strong  magnetism  as  a  consequence 
of  a  phase  transformation.  It  is  an  alloy  of  NiCoMnSn  realized  recently  in  James’ 
laboratory.  A  plot  of  magnetization  vs.  temperature  is  shown  in  Figure  4.  Note  the 
unusually  strong  magnetization  in  this  alloy,  nearly  the  strength  of  iron,  in  the  high 
temperature  austenite  phase.  A  few  degrees  temperature  change  leads  to  a  dramatic 
turning  on/off  of  strong  magnetism. 
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Fig.  6.  Magnetization  vs.  temperature  in  an  alloy  of  NiCoMnSn. 
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